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RESPONSE SOLUTIONS FOR QUASI-PERIODICALLY FORCED, 
DISSIPATIVE WAVE EQUATIONS 


RENATO C. CALLEJA, ALESSANDRA CELLETTI, LIVIA CORSI, AND RAFAEL DE LA LLAVE 

Abstract. We consider several models of nonlinear wave equations subject to very 
strong damping and quasi-periodic external forcing. This is a singular perturbation, 
since the damping is not the highest order term. We study the existence of response 
solutions (i.e., quasi-periodic solutions with the same frequency as the forcing). 

Under very general non-resonance conditions on the frequency, we show the existence 
of asymptotic expansions of the response solution; moreover, we prove that the response 
solution indeed exists and depends analytically on e (where e is the inverse of the 
coefficient multiplying the damping) for e in a complex domain, which in some cases 
includes disks tangent to the imaginary axis at the origin. In other models, we prove 
analyticity in cones of aperture tt/2 and we conjecture it is optimal. These results have 
consequences for the asymptotic expansions of the response solutions considered in the 
literature. The proof of our results relies on reformulating the problem as a fixed point 
problem, constructing an approximate solution and studying the properties of iterations 
that converge to the solutions of the fixed point problem. 
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1. Introduction 

In recent times there has been extensive interest in strongly damped systems, namely 
systems in which the term describing the damping contains a factor (where e is a 
small parameter), and subject to external forcing. Since the damping is not the term 
which corresponds to the time-derivative of highest order, this is a singular perturbation 
in £. We are interested in Ending response solutions, i.e. solutions which have the same 
frequency as the forcing term. 

A first try to understand these problems is to use perturbation theory in e and obtain 
formal series in powers of e. Nevertheless, since the perturbation is singular, one does not 



RESPONSE SOLUTIONS FOR DISSIPATIVE WAVE EQUATIONS 


3 


expect that the resulting formal series is convergent and one needs to use re-summation 
techniques to obtain that there is an analytic solution defined in an open complex domain 
which does not include e = 0, but has it on the boundary. This approach has been used 
for ode’s in |GBDn5| IGBPObl IGenlOal IGenlOb] . Different arguments for other singular 
perturbation problems can be found in [Bal94j . 

In [CCdlLlS] one can End an alternative approach for singular problems in ODE’s, 
which inspired our treatment for PDE’s. One considers the perturbative expansion to 
low orders and obtains a reasonably good approximate solution in a neighborhood of e = 0 
(i.e., an expression that solves the equation up to a small error). Then, starting from 
the approximate solution, one switches to another perturbative method (a contraction 
mapping argument) to prove the existence of a true solution. Since the problem is 
analytic in e for e ranging in a complex domain, one obtains analytic dependence in e of 
the solution for e in a certain domain which does not include any ball centered at zero. 
Indeed, we find that there are arbitrarily small values of e for which the map is not a 
contraction and the method of proof breaks down. We conjecture that this is a real effect 
and not just a shortcoming of the method. 

To motivate the procedure adopted in [CCdlLlS] . we argue heuristically that since 
£ = 0 is the most singular value of e, one attempts to do as little work as possible based 
on it. One tries to implement a perturbation theory on small but non-zero values of 
e; as soon as one gets even a flimsy foothold on non-zero values of e one switches to 
another perturbation method that is not affected by singularities (even if it contains 
some large terms, they can be beaten by pairing them with small ones). This procedure 
is somewhat reminiscent of some works in celestial mechanics, notably Hill’s theory of 
the Moon f |Hil78] . |Poi87b( Vol 2]), in which one uses a perturbation theory from an 
intermediate model which is controlled in turn by another perturbative argument. 

As it happens often in perturbative expansions, the way one deals with the hrst order 
term is different from the subsequent ones: this is even more evident in situations like 
the present one, since we are dealing with singular perturbations. In |CCdIL13] the hrst 
term of the expansion, corresponding to £ = 0, was obtained by means of an implicit 
function theorem, but the subsequent steps were all similar and they involved the same 
hypotheses. In this paper, the difference between the zeroth order term and the higher 
order ones is even more dramatic. The term in the expansion corresponding to e = 0 is 
very different from the others and in principle can be dealt with a variety of methods, 
including implicit function theorems (at least for certain cases, as we do for the model 
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described by fl2.3p below - see Section I7.2p or using variational methods (as we can do 
for the model described by (12.ip below), depending on the model we are studying. As we 
will see, when we apply variational methods, we may get even inhnitely many solutions 
of the order 0 equation. Each of them will lead to a family of solutions, which is analytic 

in e. 

Hence, the procedure adopted in the present work has two steps, with the hrst step 
having two substeps. 

a) Obtaining an approximate solution to high order, and precisely: 
al) obtaining the order zero solution; 

a2) obtaining high order approximations. 

b) Polishing off the approximate solutions to obtain true solutions. 

Each of these steps has its own methodology (indeed, step al) will be accomplished by 
means of several different methodologies depending on the model) and requires different 
conditions on the frequency as well as different non-degeneracy assumptions. Hence, the 
conditions required in the main theorem are obtained by joining together the conditions 
of all the steps. 

Nevertheless, the hnal assumptions are very weak. For example, the non-resonance 
conditions needed to carry out the whole problem are weaker than the Brjuno condition 
and they allow exponentially growing small divisors. 

The strategy above is widely applicable. In this paper we decided to document its 
breadth by applying it to 4 different models in the literature with several variations, such 
as different boundary conditions. We call these models A, A’, B, B’ (more details will 
be given in Section [2]). On the other hand, we have not optimized the hypotheses: It 
seems clear that one could obtain slightly sharper domains of analyticity, better regularity 
conditions, less assumptions on the domain, etc. We conjecture (and present arguments 
in favor) that the domains obtained are essentially optimal (see Section [H]). 

The main result for step a) is Theorem [TOl the main result for step b) is Theorem [13] 
and the hnal result is Theorem [9l 

The step al) is the solution of a functional equation. The step a2) is a Lindstedt 
procedure, which entails very mild conditions on the small divisors and requires very 
weak non-resonance conditions on the frequency. In this way one produces polynomials in 
e which solve the equation up to some (high) power of £. Under a bit stronger conditions 
on the small divisors, the Lindstedt procedure provides the existence of a formal solution 
up to all orders in e (see Theorem fTOl) . As it turns out, the solutions will be unique once 
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we fix the solution of order 0; however, as already pointed out, this solution to order 0 
can be very non-unique. 

The step b) is based on a contraction mapping principle. Hence no small divisors 
are involved but, on the other hand, we need to consider e in an appropriate complex 
domain to carry out the argument. We also note that step b) also works in cases where 
the spectrum of the operators driving the evolution is not discrete. Unfortunately, we do 
not know good conditions that ensure that one can perform step a) when the spectrum 
is not discrete. If, by any chance, one is dealing with a particular problem having a 
continuous spectrum and step a) can be performed, then step b) can be performed too 
and one can obtain the result. 

The hnal result is that the response solution is an analytic function of £ dehned in 
a domain (selected in step b) ) which does not include zero, even if it might include 
circles with real centers and tangent to the imaginary axis. Hence, the method does not 
guarantee that the Lindstedt series (the formal power expansion) converges, because the 
analyticity domain established does not contain any circle centered at the origin. Indeed, 
in |CCdlL13] one can hnd arguments that suggest that the Lindstedt series does not 
converge in general, even in the case of ODE’s. Here, we also present similar arguments 
in Section [HI 

Nevertheless, the domain of analyticity established here for models A, A’ (describing 
dissipative wave equations) is large enough, so that the application of the Nevanlinna- 
Sokal theory f [Nevl9[ ISok80[ IHar49] ) on asymptotic expansions applies. As a conse¬ 
quence, the response solutions constructed here have an asymptotic expansion and these 
functions can be reconstructed from their asymptotic expansions by re-summation. No¬ 
tice that this procedure is very different from establishing the existence of the solution by 
re-summing the series. Of course, since the problem is nonlinear, re-summing the series 
is not enough and one needs other arguments to show that the re-summation solves the 


equation ( |Har49j , see also |Bal94[ IBLSn2[ IGBDn5[ IGBPObf IGenlOaf IGenlOb] ). 

In some models such as models B, B’ (describing large stiffness equations), we obtain 
domains of analyticity which are cones containing the imaginary axis and have an aper¬ 
ture of 7r/2. We conjecture that these domains are essentially optimal (see Section [H]). 
We will show that the functions we construct have the same asymptotic expansions as the 
formal power series. On the other hand we note that in domains of this kind it is not clear 
that the response solution can be obtained by re-summing the asymptotic expansions: 
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indeed in these domains there are non-trivial functions whose asymptotic expansion van¬ 
ishes, so that the expansion is not unique (e.g., the Cauchy example exp(—which has 
an asymptotic expansion vanishing in domains of aperture vr/2). As a consequence, it 
could well happen that for these models the response solutions lead to exponentially small 
phenomena. Notice that model B is an inhnite dimensional analogue of fast oscillators 
for which exponentially small phenomena have been established (see [BFGSl^ i. 

In |CCdlL13] the problem considered is the varactor equation, which is a single ODE. 
Even if the step b) in |CCdlL13] was just an elementary one (based on contraction argu¬ 
ments), the results obtained in |CCdlLl,^ improved the existence domains and weakened 
the non-resonance conditions that have been imposed in the previous literature. It seems 
plausible that using a more efficient hxed point argument (e.g., a KAM theory) or higher 
order perturbations in step a) would improve the results. The analyticity domain has 
later been extended for ODE’s in |CFG13] (where a domain of analyticity tangent more 
than quadratically to the origin was established) and the non-degeneracy assumption 
on the non-linearity has been relaxed (for real e) in |GFG14] . As further references, we 
mention also |Rab67l IRab68( IGra83] . where the periodic case with real small damping 
has been considered. 


1.1. Description of the main results. The goal of this paper is to extend the method 
of |GGdlLL^ to some PDE’s. The method is very flexible and we will present results 
for four different models considered in the literature, each with three different types 
of boundary conditions (see Section [2]). It is clear that there are many more models 
that could have been considered by the method. Of course the main difficulty of the 
extension to PDE’s is that the operators are unbounded. Hence the reformulation of 
the problem as a fixed point problem requires some more thought, even to get a viable 
formulation. For example, we need to ensure that the operator maps some space into 
itself and that the space satisfies suitable properties (such as Banach algebra properties). 
So, considerable effort goes into the choice of spaces as it happened in the classical study 
of elliptic problems (see Section [33]) . 

The models we consider in this paper have the form 

duu{t, x) H—Friction — A^u{t, x) -|- h{u{t, x),x) = f{ojt, x) 

£ 

for models A, A’ below, or 

x) -|- dtu{t, x) — A^uit, x) + non-linearity = /(cat, x) , 
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for models B, B’ below. The equations will be supplemented with the boundary condi¬ 
tions. 

In all models, given a domain V (as specihed in Section 12.ip and denoting by V its 
topological closure, n : M x P —?■ M is the unknown. We will require that the following 
data of the problem are hxed: 

• The boundary conditions; 

• /i:MxD—)"Mto which we refer as the non-linearity; 

• / : X D —)■ R (with T'^ = (R/Z)'^) to which we refer as the forcing; 

• cj G R'^, which denotes the frequency of the forcing. We assume without loss of 
generality that ca has rationally independent components, namely that: ca • /c 7 ^ 0 
for all A: G Z'^ \ {0}. 

Of course we assume that the forcing and the non-linearity are such that the boundary 
conditions are maintained. We will also need: 

• Quantitative estimates on the size of \u: ■ k\~^ as a function of \k\ (which will turn 
out to be weaker than the Diophantine or Bryuno conditions); 

• A non-degeneracy condition on the non-linearity. 

Then we shall prove the following “meta”-result. 

Meta-Theorem 1. Under the above requirements there exists a response solution for 
the models of the form above, defined for e in an appropriate complex domain: the spe¬ 
cific form of such domain depends on the model considered as well as on the boundary 
conditions. 

The precise statement of the result requires the introduction of the spaces, the domains 
and a precise formulation of the regularity condition that we will give later on, see 
Sections [SHU The existence of the solutions of the zeroth order term is discussed in 
Section[5]and Appendix[Bl The proof for the case of dissipative wave equations is provided 
in Section [6], while the modihcations of the proof for the other models are given in 
Section [71 Some arguments supporting that the domains are almost optimal are given in 
Section [HI 


2. Models considered and some preliminary assumptions 

In this section we present the models we intend to study and we state the required 
non-degeneracy assumptions. 
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In what follows we will assume that A^; is a self-adjoint elliptic operator of second 
order; in the physical applications we have in mind it is the Laplace-Beltrami operator. 
We will not necessarily assume that A^. is a constant coefficient operator. 


2.1. PDE’s considered. We will consider PDE’s for which the space variables range 
in the topological closure of a domain P and we will look for solutions quasi-periodic in 
time. The domain P can be: 

Dl) a compact manifold without boundary, for example (we will refer to this as 
the periodic case), 

D2) an open, bounded, connected subset of with a C°° boundary. In this case, we 
will supplement the solutions with either Dirichlet or Neumann boundary conditions. 

Therefore, we will consider the following (standard) boundary conditions, each one 
leading to a different functional setting, which we will specify below: 

D) Dirichlet boundary conditions, 

N) Nenmann boundary conditions, 

P) Periodic bonndary conditions. 

Following the usual practice, we interpret the boundary conditions as describing a space 
of solutions: the operator A^. acts on this space and of course the spectral properties of 
Aa; depend on the space too. Of course, in order to specify the function space we also 
need to specify a norm. Our treatment will be for spaces of functions which are analytic 
in t and differentiable in x. 

We will consider four different PDE’s: models A, A’, B, B’ below. Each of them may 
have entirely different boundary conditions (Periodic, Dirichlet and Neumann). 

A) The dissipative wave model: The hrst model is a direct analogue of the 
varactor equation studied, e.g., in [CCdlL13l ICGV05( IGBD05( IGBD06( IGenlOa) IGenlOb( 
IGFG141 IGFG13j : the model is obtained from the wave equation by adding a singular 
friction proportional to the velocity: 

dttu(t, x) + -dtu{t, x) — A^u{t, x) + h{u{t, x),x) = f{ujt, x) . 


( 2 . 1 ) 
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A’) The frequency over-damped model: We modify the friction of model A, as 
described by the following eqnation 


1 


dttu{t, x) H— dtAxu{t, x) — Axu{t, x) + h{u{t, x),x) = fiut, x) 


( 2 , 2 ) 


In this model, which has been studied for instance in |PSMn9] . the damping is stronger 
for the spatial modes with larger spatial frequency. Indeed, not only the damping term 
dfA^u in fl 2 . 2 l) is affected by a factor which is the inverse of the small parameter e, 
but it contains also the unbounded operator A^. 

For simplicity, we have considered the case where the A^. appearing in the damping 
and in the restoring force are the same operator. Some slight generalizations are possible, 
such as taking different operators for the damping and the restoring force provided they 
commute. In many physical applications, it is natural that the operators describing the 
damping and the force commute, since they have to be translation invariant and isotropic. 

B) Large stiffness model: This is a generalization of the model introduced in 
[FMillFMPSnr , described by the equation 


e‘^duu{t, x) + dtu{t, x) — A^uit, x) + h{u{t, x),x) = f{u!t, x) ; 


(2.3) 


in [Flol4[[FMPSOTj one can hnd the specihc case h{u, x) = 7 /( 1 with e G M and with 
7 > 0 a dimensionless parameter which provides the relative strengths of electrostatic 
and mechanical forces. 

Equation (12.31) models an electrostatically actuated MEMS (Micro-Electro-Mechanical- 
Systems) device. Precisely, the physical interpretation of the model (12.31) is that the 
restoring force of the oscillators forming the wave equations is very large. This type of 
equations are used to model the deflection of an elastic membrane suspended above a 
rigid ground plate, with a voltage source and a hxed capacitor. The model represents 
the limit of small aspect ratio, when the gap size is small compared to the device length. 
The paper [FMPSOT] contains a detailed discussion of the motivation. It is interesting to 
note that the varactor equation is somehow a model for the problems considered here. 

B’) The modified large stiffness model: We modify the non-linearity of model B 
by assuming that it is of order e, as described by the following equation; 


e'^duuit, x) -|- dtu{t, x) — A^uit, x) + eh{u(t, x),x) = x) ; 


the above equation appears in the study of MEMS with high aspect-ratios and/or when 
the applied tension is high. 
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2.2. Regularity assumptions and boundary conditions. We will require that / is 
smooth in x, satishes the boundary conditions, and is analytic in the variable 6 = ut. We 
will formulate this assumption more precisely by saying that / belongs to a Hilbert space 
which we shall call Apj^m] see the definition in Section [3.5.21 where we will impose some 
restrictions on the parameters (as we will see p, j measure the analyticity properties and 
m measures the regularity properties in the space variables). 

We will assume that h has some regularity properties too. Roughly, we will require that 
h is analytic in its first argument and differentiable when the second argument ranges over 
V. Slightly more precisely, we will require that h is such that given a function u G Ap^j^rn, 
then h{u{9,x),x) is also in Apj^m and that the map u i—)■ •) is differentiable in the 

sense of maps in Banach spaces. We will also require that h satishes certain geometric 
conditions ensuring that the boundary conditions are preserved. Precisely we make the 
following requirements. 

BCD. For Dirichlet boundary conditions we require that h{0,x) = 0. 

BCN. For Neumann boundary conditions we require that 

n{x) ■ {Dxh){u,x) = 0 for all x G dV , u G M , (2.4) 

where n{x) denotes the normal to the domain V at x. In this way, we obtain that 

n{x) ■ Dx[h{u{t, x), x)] = {Duh){u{t, x), x) n{x) ■ D^uit, x) + n(x) • {Dxh){u{t, x), x) , 

which equals zero if u(t, •) satishes the Neumann boundary conditions and fl2.4l) holds. 
Notice that we have used that D^h is one-dimensional. 

We anticipate that, besides the above regularity and boundary conditions, we will also 
require some non-degeneracy conditions on h. 

Remark 2. In this paper we will construct solutions analytic in time. The proofs work 
similarly in spaces of functions with Sobolev regularity in time (with high enough Sobolev 
exponent depending on the dimension of the freguency), when developing the theory for 
finitely differentiable cases (i.e., when the functions f,h are only assumed to be finitely 
differentiable) . Of course in this case one can consider only e G M. 

3. Formulation of the problem and overview of the method for model A 

In this section we go over the method for model A and reduce it to a hxed point 
problem. We will present hrst the formal manipulations, since they are the motivation 
for the constructions and the precise dehnitions given later on. Notably, the choice of 
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spaces in Section 13.51 will be motivated by the need that the operator appearing in the 
hxed point equation maps the spaces into themselves and it is a contraction. 


3.1. Response solutions and formal power series. Our goal is to hnd response 
solutions of the form 


Ueit,x) = Co{x) + Ue{ujt,x) , 


(3.1) 


where for each hxed e, t4 • x P — )■ M is at least 0{e). We will refer to cq as the zeroth 
order term and omit the index e whenever this does not lead to confusion. 

We will hrst show that when we write 17 as a formal power series in e: 


u, = Y.^>u, , 

i=i 


(3,2) 


the coefficients Uj can be formally dehned: the appropriate Banach spaces of functions in 
which the coefficients actually exist will be specihed in Section 13.51 Such Banach spaces 
will include regularity properties as well as the boundary conditions. 

Inserting (13.Ih in (12.11) . we get that the function must satisfy the equatior0: 


(3.3) 

The solution of (13. 3 h will be the centerpiece of our treatment. Later, we will develop 
analogous procedures for models A’, B, B’ (see Section [7]). We remark that the series 
expansion fl3.2p does not contain the term j = 0; in fact, if we add a term Uq to the series 
(13.2p . then taking the coefficient of order in (I3.3p . the term Uq would satisfy 


(cj • S/b)Uq{0,x) = 0 , 


showing that the solution t/o, which can be found under the non-resonance assumption 
on cj, is independent on 9. However, having written the response function as in fl3.ip 
with Co being the 6*-independent part, we conclude that it must be Uq = 0. 


3.2. Formal solutions of the equation for response functions. In this section we 
describe how to obtain a formal power series solution for fl3.3p . This is step a) of the 
strategy discussed in the introduction. 

^The search of quasi-periodic solutions with frequency w having rationally independent components 
is equivalent to looking for a solution u = u{6,x) of the differential equation in which ixt is replaced by 
9 and dt is replaced by x ■ Vg; this is why we shall study functions of the form u = u{6, x). 
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3.2.1. Dividing the problem into zeroth order and higher orders. We introduce the nota¬ 
tion: 

N,U{e,x) = [{cu-Vgf + ^{uj-Ve) + C]U{e,x), (3.4) 

Cr]{x) = -A^r]{x) + h'{co{x),x)r]{x) , (3.5) 

G{U){9,x) = h{cQ{x) + U{9, x), x) — h{cQ{x), x) — h\cQ{x), x) U{9, x) . (3.6) 

Note that the operator depends on e, whereas G and C are independent of £. 

Remark 3. If the operator C is elliptic and self-adjoint in L\c{T^) (namely Lf{D) with 
boundary conditions), then the eigenfunctions constitute a complete set for the Hilbert 
space L|p(P) because C has compact resolvent, see |Hel69( [Dav95j . Similar considera¬ 
tions apply to the analogous operator C introduced for models A B. 

We will use that the eigenvalues of C are real and that we can characterize the 
Sobolev spaces in terms of the coefficients of the eigenfunction expansions. 

With the above notations and denoting by (•) the average with respect to 9, it is just 
elementary algebra to show that the equation fl3.3p is implied by the pair of equations: 

NM9,x) + G{U,){9,x) = f{9,x) - (/)(x) , (3.7) 

- Aa:Co(x) h{co{x),x) = {f){x) . (3.8) 

The reason to divide the equation fl3.3p into fl3.7p and (13.81) is that (13.8p is the leading 
order in e. 

Notice that the order e'^-term in (13.31) is 

{u ■ Vo)Ui{9, x) - A,,co{x) + h{co{x),x) = f{9, x) . (3.9) 

Hence for a solution Ui of (13.9p to exist, it is necessary that the average of (13.9p with 
respect to 9 is zero (hence equation (I3.8p j. Of course, if u satishes suitable non-resonance 
conditions and the functions are smooth, it is indeed possible to obtain Ui. The solution 
of equation (13.9p . which is standard in KAM theory and which can be dealt with Fourier 
expansions, will be discussed in Section 16.11 In conclusion the system (13.7p . (13.8p is 
equivalent to (13.3p . if we look for formal solutions as in (13.21) . 

Notice that the system (13.71) . (13.8p has an upper triangular structure. In particular, the 
equation (13.8p involves only cq: once we obtain a solution cq of (13.8p . we can substitute 
it in (13.91) and obtain the solution Ui, and then proceed to higher orders. 
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The existence of solutions of fl3.8p has been studied extensively in the literature through 
a great variety of methods. In Appendix [B] we will present some of the results available 
in the literature. 

As a result we obtain that under many circumstances there are several (often inhnitely 
many) cq solving (13.Sh . For each of them we will see that (under appropriate non¬ 
degeneracy conditions) we can hnd a unique solution (hrst as formal power series and 
then as analytic function in a domain). Hence the upper triangular system may have 
many solutions, but the only source of non-uniqueness is the equation (13.Sh for cq. 

3.2.2. Preliminary assumptions on the operator C. We now specify the spectral proper¬ 
ties of the operator C. In general, to characterize the spectrum of an operator, one needs 
to specify on which space it acts. Nevertheless in our case we assume that the operator 
is elliptic and that the domain is compact. In this case the spectrum is discrete and it is 
the same in all Sobolev spaces. 

The well known reason why the spectrum is independent on the spaces ( |Kat76l IHel69] ) 
is that, when the operator is elliptic, the solutions gain regularity: this translates to the 
fact that the resolvent is compact on any Sobolev space and hence, for all the Sobolev 
spaces the spectrum is just a discrete set of eigenvalues with hnite multiplicity. Further¬ 
more, again by regularity theory, the eigenfunctions are very smooth, so that they are 
eigenfunctions in all Sobolev spaces. Therefore the spectrum is the same in all Sobolev 
spaces. The following assumptions H1-H2 will be requested for model A as well as for 
models A’, B, B’ once the operator C is suitably dehned. 

HI The spectrum of C is discrete and its eigenvalues satisfy: 

0 < A„ < A„+i , V n > 1 

and the multiplicity of each eigenvalue is hnite, possibly increasing with n; 

H2 The smallest eigenvalue is positive: Ai > 0. 

In the case of models A’, B’ we will also assume the following hypotheses on the 
operator —A^.. 

HI’ The spectrum of — A^, is discrete and its eigenvalues A^ satisfy 

0<A^<A^+i, Vn>l 

and the multiplicity of each eigenvalue is hnite, possibly increasing with n. 

H2’ The smallest eigenvalue is positive: Af > 0. 
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Remark 4. Note that a consequence of HI and HI’ is that there is an orthonormal 
basis of eigenfunctions for the operator P = C or — in L‘^{V), such that 

P^n = for n = 1 , 2 ,with Af ^ = 

Remark 5. We note that HI and H2 are in turn assumptions on h, cq. In some 
arguments, we will need to assume only HI, but in order to get the crucial estimates on 
the “small divisors” (and hence to obtain the final result, see Theorem lTSi) . we will need 
to assume that there is the spectral gap in H2. 

The above assumptions can be slightly modihed; in particular, the previous assump¬ 
tions H1-H2-H1’-H2’ can be extended to encompass the case of a continuous spectrum 
(see Remark [32] below). 


An if P = £ , 

A^ if P = -A, . 


3.2.3. The nonlinear term and the boundary conditions. We start by noticing that G is 
the functional analogue of the nonlinear term used in [CCdlL13] . provided of course that 
the operator G is dehned from some appropriate space to itself. 

In this section we will just check that, if we assume h to satisfy the conditions BCD, 
BCN in Section 12.21 then the operator G preserves the spaces of functions satisfying 
these conditions. 

In the case of periodic boundary conditions, there is nothing to check. 

For Dirichlet boundary conditions we observe that if x G dP and cq, U satisfy the 
Dirichlet boundary conditions, then Co{x) = 0, U{6,x) = 0 and, hence G{U){6,x) = 0. 

For Neumann boundary conditions we observe that, if h satishes fl2.4l) . then we just 
need to check that 


n{x) ■ D^[h'{co{x),x)U{e, x)] =0 Vx G 0 G . (3.10) 

The left hand side of (13.101) can be written as the sum of three pieces, i.e. for x G dT>, 
n{x)-D„[h'{co{x),x)U{9,x)] = h''{co{x),x){n{x) • P^co(x))P(6 ',x) 


n{x) ■ {D,,h'){co{x),x)U{9,x) + h'{co{x),x) n{x) ■ {DJJ){9,x) . 


(3.11) 


The hrst and second terms in the right hand side of (13.111) vanish, since we impose that 
Co(-) satishes Neumann boundary conditions and h satishes (12.41) . Therefore, if U{9,-) 
satishes Neumann boundary conditions, the last term in (13.lip will also be equal to zero. 
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3.3. The higher order equations. We are looking for a formal power series solution 
of (I3.7jl . Assume that we solved (13.Sp . insert (13.2p into (13.7p and expand the power series 
(this requires enough regularity for the function h which we will make explicit later). 

Equating the coefficients of the same power for > 0, we obtain the following 
recursive equations for > 0: 

(u;-Ve)t/^+i( 0 ,a:) + {u-VgfUN{9,x) - ^^Un{9,x) + h'{c,{x),x) Un{9,x) 

= ^^(co(0,a;),17i(0,x),...,17^_i(0,x)) , (3.12) 

where Sn is a polynomial expression in 17i,..., Un-i obtained by applying the Taylor 
theorem to order N in the equation (13.3p and gathering terms. 

We think of (13.121) as an equation for f/Ar+i, given all the previous terms of the expan¬ 
sion. Of course we need to assume that ca • /c 7 ^ 0 and indeed that it is not too small as 
\k\ increases. Provided that 

{A^Un{9, x) - h'{co{x),x) Un{9, x) + S'Ar(co, Ui,..., Un-i){9, x)) = 0 , (3.13) 

we can hud I/at+i which is unique up to the choice of an additive function of x alone. 

Hence, as it is standard when dealing with Lindstedt series, proceeding by induction 
we assume that we have determined f/i,... Un and then using (I3.12p we can determine 
T/tv+i up to an additive function of x\ such a function is obtained by solving (I3.13P and 
this can be done because of H2. 

Sufficient conditions for the existence of an approximate solution provided by a trunca¬ 
tion to order N of the series expansion (13.21) are given in Theorem 16.1( see also Section iQ] 
for a discussion of the existence of an approximate solution to a finite order by solving 
(I3.12P (compare with (16.21) in Section [ 6 H]). 

3.4. Formulation of the fixed point problem equivalent to (13.31) . As we shall see, 
the operator in (13.41) is invertible in the spaces Apj^m alluded above, if £ ranges in a 
suitable domain, so that (13.7p can be rewritten as 

U,{9,x)=N-^ [-G{U,){9,x) + f{9,x)-{f){x)]^%{U,)i9,x) , (3.14) 

where we have introduced for convenience the operator %] we will show that (13.141) can 
be solved by a contraction mapping argument. 

Therefore one of the crucial points of the strategy will be to study the invertibility of 
and give quantitative estimates on its inverse, notably the Lipschitz constants. In 
order to do so, we provide a uniform lower bound on the eigenvalues of (which will 
depend on e), using the assumption H2 on the eigenvalues of C. By carefully examining 
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such ^-dependent bounds, we will show that, for e in a suitable domain, the operator 
appearing in the right hand side of (13.1411 sends a ball centered at the approximate 
solution (given by the perturbative expansion) into itself and that it is a contraction 
inside this ball. Hence, the hxed point can be obtained by iteration, starting from the 
approximate solutions. 

We think at the iterative procedure as taking a function analytic in e and producing 
another analytic function of e. We will show that the convergence is uniform for e in a 
suitably chosen complex domain. Then it is a standard argument that the limit is an 
analytic function of e in this domain. 

The contraction mapping argument is classical; however it requires to use spaces in 
which we have sharp estimates, so that we do not lose any regularity and we obtain that 
the operator in (13.1411 sends the spaces into themselves. 

Of course, once we have dehned the spaces, we will have to justify the formal manipula¬ 
tions, such as the existence of functional derivatives. This amounts to making regularity 
assumptions on the term h, which justify the use of the Taylor’s theorem up to order N 
for the composition operator. 


3.5. Choice of spaces. In this section we present the spaces we will use. We discuss 
some of their elementary properties in Appendix where we also add a remark about 
the continuous spectrum. 

The leading principle is that the norms of the functions can be expressed in terms 
of generalized Fourier coefficients, namely the coefficients associated to the basis given 
by the product of the Fourier basis in 6 and the eigenfunctions of C with boundary 
conditions in x. 

This principle allows us to estimate rather easily the inverse of the linear operator 
W in (!3.4jl just by estimating its eigenvalues, because we are allowed to use the base in 
which Ag is diagonal. 

We also need the spaces to have other properties allowing us to control the non-linear 
terms, such as Banach algebra properties and properties of the composition operator, so 
that we can study the operator G. Since we want to obtain analyticity in e, we will also 
need spaces of analytic functions and, in order to simplify the analysis, we require that 
they are Hilbert spaces. Note that we think of the functions in x as “scalars” in analogy 
to what happens in [CCdlLl,^ : hence, it is natural to consider Hilbert spaces of analytic 
functions in 6 taking values in another Hilbert space of functions of x. 
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The choice of the spaces presented here satishes such properties and leads to simple 
proofs. Of course we are not claiming that the choices we make are optimal and it is 
quite plausible that other choices (e.g., analytic functions in both variables) could lead to 
better regularity. The main problem in using spaces of analytic functions in x is that it 
is not clear to us how to express the analyticity of a function in terms of the coefficients 
of the expansions in eigenvalues. 

We will present several equivalent norms, since some of the properties of the space will 
be easier to verify in one norm than in another. Henceforth, given two (hnite or inhnite 
dimensional) equivalent norms || • ||, |1 ■ |1^ we write || • || = || • ||^ . 

3.5.1. Sobolev spaces with boundary conditions. In this section we introduce the Sobolev- 
like spaces which we will use; we will dehne them only for indices m G 2N, since this is 
enough for our purposes. The advantage is that, for these indices, it is possible to give 
particularly simple characterizations of the norm in terms of the eigenfunction expan¬ 
sions. Using several characterizations of the norms allows one to obtain simple proofs of 
Lipschitz properties of operators. 

For functions S' : P —)■ C satisfying the corresponding boundary conditions and for 
m G 2N, we dehne the family of equivalent norms as 

. (3.15) 

If S{x) = Sn^n{x) with S'n G M and as in Remark 01 then the Sobolev norm 

(I3.15P is given by 

OO 

iisiii,j. = 5^A:ris„iu 

n=l 

where A„ are the eigenvalues of C (recall that the $n’s form a basis of eigenfunctions of 
C). Since C is elliptic, by Garding’s inequality (see |Taylla| -Theorem 6.1 of Chap. 7), 
we have 

ll'S'||i72‘ “ II'^IIr’" ’ (3.16) 

where || • \\h^ is the standard Sobolev norm, namely 

||F||^,^ = ||(Ao + 1 )™/2R|U2 

with Aq the standard constant coefficient Laplacian and we are considering S satisfying 
the specihed boundary conditions. 

The spaces H^{V)bc and H'^{V)bc are the completion of - the set of C°° func¬ 
tions with compact support contained in the interior of "D - under the above norms. 
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For notational convenience we shall not write explicitly the dependence on the bonndary 
conditions nnless needed. 

It is well known that for m > £/2 the Sobolev spaces satisfy the Banach algebra 
property ( [Tayllb] ) and hence the equivalence of the norms || • \\h^ and || ■ \\h^ in (Id.lbh 
implies for every Si, S 2 G 

for some constant C > 0. 

When m > i/2 the Sobolev embedding theorem says that the functions in H'^ are 
continuous, so that the Dirichlet boundary conditions have classical meaning. 

Similarly, when m > i/2 + 1, the gradient of functions in are continuously differ¬ 
entiable. Hence, the Neumann boundary conditions have classical meaning. 

3.5.2. Spaces of analytic functions of complex variables taking values into Banach spaces. 
We introduce domains that consist of a strip around the torus in the imaginary 
direction. We will consider analytic functions in these domains. 

Definition 6. Given p > 0, we denote by the set 

= {0 e (C/Z)" : Re(0,) e T , |lm(0,)|<p, j = 

When we consider functions of 6^ G and x G P, we can think of them as functions 
from into which are analytic §. The spaces which we will consider are the standard 
Bargmann spaces taking values into 

Given a function u = u{6,x) which we expand as 

u{9,x) = ^ e‘^^"’^'^^n(,x)uk,n , (3.17) 

we will consider the space of analytic functions of 6 endowed with the H"^) norm 

dehned below. We emphasize that we are considering as a 2d-dimensional real manifold 
with boundary. Again, for simplicity, we just consider the even Sobolev exponents j. 

^ When we consider domains which are closed with a smooth boundary, we refer to analytic functions 
as functions which are analytic in the interior and that extend continuously to the boundary. For us, 
using domains which are compact is slightly more convenient in order to quote embedding theorems, 
etc. Nevertheless, in order to avoid repetitions, we omit that analyticity is meant only for the interior 
and that we assume the extension to the boundary. 
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Precisely, for p > 0, j,m E 2N, setting (where the bar denotes 

complex conjugation), we dehne the norm as: 



(3.18) 


We denote by Apj^m the space of functions analytic in 6 whose norm || • ||p,ym is hnite. 
Note that Apj^m are Hilbert spaces, since the norm flS.lSp clearly comes from the inner 
product 


{u,v) = / (m, (Ag + iyv)Hrr^ d‘^‘^9 . 

Jti 

Moreover they are complete, since the limit in the || ■ ||pj,m-norm of analytic functions is 
an analytic function f |RS80j h 

Remark 7. We think of a function u G Apj^m o^s an analytic function from into if™, 
say 6 -E u{9, •). In this way, the problems considered here look closer to the formulation 
of the varactor problem considered in [CCdlLlS] , The PDE looks formally like an ODE 
in if™ and the response solutions will be analytic functions from the torus into hi™. Eor 
sufficiently high m, these will be classical functions which are analytic in the t variable 
and differentiable in the variable x. Hence, they will be classical solutions for the PDE. 


4. Precise statement of the results 


4.1. Approximate solutions of the fixed point problem. For all models we can 
give a dehnition of “approximate solution” as follows. 

Definition 8. Let us consider a family of functional equations 




(4.1) 


where He '■ Apj^m —t Apj^m is an operator that maps e-dependent families into families 
(of course the operator He may have an explicit e-dependence). We say that 


M 



for some M G is an approximate solution up to order M of fl4.1l) . if 


(4.2) 
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4.2. Main results. Our main results are provided by the following Theorems [TOl [T3l 

Theorem [9] is based on a contraction mapping argument and it states the existence of 
a solution, provided we assume the existence of an approximate solution. 

Theorem [10] gives sufficient conditions for the existence of an approximate solution up 
to any order under some non-resonance assumptions on the frequency: the higher is the 
order of approximation we want to reach, the more restrictive will be the condition on 
the frequency. 

Theorem [13] summarizes the results above, i.e. it gives the existence of an analytic 
solution under the requirements of Theorem [TO] which provides the approximate solution. 
The proof of Theorem [13] relies on applying Theorem ]9] to the approximate solutions 
provided by Theorem [10] 

In the Theorems [9] [IHl [I3] we will assume that / belongs to the space of functions 
Apj^m as in Proposition ]3ni which ensures the validity of the Banach algebra property. 

Theorem 9. Assume that f is in Apj^m for p > 0, j,m E 2N, j > d, m > ^|2 
(m > Ij^ in the case of Neumann boundary conditions). Let h : B x V ^ C with 
B G C open set, and let V he either of the form D1 or D2 as in Section\2fJ\ We assume 
that h is analytic in B and C^iV) fl C{V) in x. 

Consider the models A, A’, B, B’ with D, N, P boundary conditions; assume that the 
hypotheses H1-H2 are satisfied (see Section VS.2.S\) . and that the non-linearity h satisfies 
BCD or BCN (See Section \2.^} in case of D or N boundary conditions, respectively 
(depending on the boundary conditions considered for the equation). For models A’, B’ 
assume also Hl’-H2’. 

For model A assume that the zeroth order term cq (see fl3.8p ) admits a solution (some 
sufficient conditions are given in Appendixand that for some M G M, M > 2, there 
exists an approximate solution in e of (13.3h up to order M. 

Let e he in the domain LIb = Uo-Ho-.b with B > Bq for some Bq > 0 sufficiently large, 
cr > 0 sufficiently small, where 

Qcr,B = {e = ^ ip E C : ^ > B , a < \e\ < 2a} (4.3) 

and 6 in the strip of size p > 0 

= {0 e {C/Zy : Re(0,) G T , |lm(0,)| <p, j = l,...,d} . 

Then, there exists a function = Lf^{6,x) E Apj^m, which provides an exact solution of 
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For model A ’ assume that the zeroth order term cq (see (13.Sh ) admits a solution (some 
sufficient conditions are given in Appendixand that for some M G N, M > 2, there 
exists an approximate solution in e of 

{uj-V 0 fUe{ 6 ,x)+-{uj-Ve) A^Ue{e,x)-A^Ue{ 6 ,x)-A^co{x)+h{co{x)+Ue{ 6 ,x),x) = f{ 6 ,x) . 

(4.4) 

Let e be in a domain of the form fl4.3p . Then, there exists a function u = u{ 6 ,x) = 
Co(x) + Us{0,x) as in fl3.ip . belonging to which provides an exact solution of fl4.4p . 

For model B we assume that there exists an approximate solution of 

s\u ■ VefUeie, x) + {cu- Vg)Ue{e, x) - x) + /i(f/,(0, x),x) = /(0, x) (4.5) 

up to order M with M G N, M > 2. Assuming that e belongs to the domain 

ns = {e = f + ir]eC : Re(-e2) > 5} U {e = ^ G M : 5 < |^| < 25} (4.6) 

for some 5 > 0, then there exists a function Us = Us{9,x) G Apj^m, which provides an 
exact solution of fl4.5p . 

For model B’ assume that the zeroth order term admits a solution (see Section \5.A} 
and that for M G N, M > 2, there exists an approximate solution up to order M in 
Apj^m of 

e^(u ■ VgfUs(9, x) + (u ■ Vg)Us(9, x) - A,,Us{9, x) + eh{Us{ 6 , x),x) = f{9, x) . (4.7) 

Assuming that e belongs to Fts as in fl4.6p for some 5 > 0, then there exists a function 
Us = Us{9,x) G Apj^m, which provides an exact solution of fl4.7p . 

In all the cases above, the solution Us is analytic in the considered domains as a 
function of e and it is asymptotic to the approximate solution. 

Note that Theorem [H] involves mainly regularity assumptions and the requirement that 
there exist approximate solutions at least of order 2. Sufficient conditions for the existence 
of an approximate solution (given by the expansion to any arbitrary order) are provided 
by the following result. 

Theorem 10. Assume that f is in Apj^m for p > 0, j,m E 2N, j > d, m > 1(2 
(m > 1(2 + \ in the case of Neumann boundary conditions). Let h : B x F ^ C with 
B G C open set, and let B be either of the form D1 or D2. We assume that h is analytic 
in B and C^{V) n C(V) in x. 

Consider the models A, A’, B’ with either D, N, P boundary conditions and assume that 
h satisfies either BCD, BCN in case of D, N boundary conditions, respectively (depending 
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on the boundary condition considered for the equation). Assume that the zeroth order 
term admits a solution (see, respectively, Appendix\M for models A, A’ and Section COI 
for model B’). Furthermore: 

i) Assume that there exists M G N, such that 

|i|-‘log|wfc|-‘<^ VfceZAD}. (4.8) 

Then, there exists an approximate solution in e of fl3.3p . (14.41) . (14.7p up to order M. In 
particular, if 

limsup logdo; •= 0 , (4.9) 

\k\^oo 

then we can obtain a formal power series in e (whose coefficients are well defined) solving 
the equation up to all orders. 

ii) If we assume that f is a trigonometric polynomial, then there is a formal power 
series in e which is a solution of (13.31) . (14.4p . (14.7p up to all orders, without requiring 
any non-resonance bound on the frequency oj. 

Remark 11. Recall that we are assuming that uj ■ k = t), k & lA, implies fc = 0. 

Remark 12. It seems likely that the condition (14.81) (which is even weaker than the 
Bryuno condition) is optimal. Each step of the computation of the perturbative expansion 
involves solving a differential equation with {u ■ k)~^ as small divisors. Hence, we expect 
that the solutions, in general, lose a domain of definition of size p/M at each step. 

Combining Theorems 191 and fTOl we obtain the following result. 

Theorem 13. Assume that f, h satisfy the assumptions of Theorem [H regarding the 
regularity, the boundary conditions and the existence of the zeroth order solution. Fix 
M >2 and assume (14.81) . Then, we have the following results. 

For model A there exists a solution of equation (13.3p . which is analytic in e and 9, and 
satisfies the D, N or P boundary conditions. The analytic solution exists for e in the 
domain Qb as in Theorem [H and 6 inT^. 

For model A’ assuming HI’, H2’, there exists an analytic solution of (14.41) m 0 G 
and e with e in a domain of the form Qb as in Theorem\B 

For model B with D, N, P boundary conditions, provided the existence of an approxi¬ 
mate solution, there exists an analytic solution of (14.51) in 9 E and e in a domain of 
the form as in (14.6p for some 5 > 0. 
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For model B’ with D, N, P boundary conditions, assuming there exists an 

analytic solution of (14.7p in 9 E and e in a domain of the form fig as in (I4.6jl for 
some (5 > 0. 


Remark 14. Note that the non-resonance condition which we need to impose on u to 
obtain the existence to all orders is more restrictive than the non-resonance condition we 
need to obtain the existence of an analytic solution defined in the domain Qa,B or 
Since, as we argued in Remark\T^ we believe that the conditions are optimal, it seems 
that given an Mq and an uj that satisfy fl4.8p for M = Mq, but not for M = Mq + 1, then 
we can obtain functions that have an expansion up to order Mq, but not Mq + 1. 

It seems, therefore, possible to arrange the existence of models with a solution analytic 
in a domain of the form orfls, but without Taylor expansion beyond a certain order. 

Remark 15. By restricting the domain in Theorem UR we can obtain stronger contrac¬ 
tion properties for the operator T defined in (I3.14p . 

For example, for model A one possibility would be to consider only one of the do¬ 
mains defined in fl4.3p with a small enough. Another possibility is to consider conic 
domains defined as 

'^ 6,(7 = {s E C : I Im(£)|/|e| < 5 , a < |e:| < 2a} 

for some 6, a > 0. We refer to |CCdlL13] for further details on the study of the solution 
of a forced strongly dissipative ODE on conic domains. By restricting to the real line, it 
seems possible to obtain results for finitely differentiable non-linearities. 


5. Existence of solutions of the zeroth order term 


The existence of solutions of the zeroth order equation for models A, A’ has a very 
extensive literature and can be done by a variety of methods; with reference to classical 
textbooks like [StrOSP lAMOTP IPrel3] , we defer the presentation of some of such results in 
Appendix [HI For models A, A’, we will produce several (even infinitely many) solutions 
of the zeroth order equation. 

In this section we confine ourselves to discussing models B, B’, using arguments based 
on an implicit function theorem in Banach spaces. However, we mention that there are 
other possibilities which we have not covered, for example methods based on index theory 
( |FG951 IBer73j ). We will show that each of these solutions of the zeroth order equation 
continues into a formal solution to all orders and that, furthermore, it can be modified 
to be a true solution. 
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5.1. Solution of the zeroth order term for model B. The zeroth order equation 
for model B is: 


(cj • S/e)Uo{e, x) - x) + h{Uo{9, x),x) = f{9, x) . (5.1) 

We will solve equation (15.Ih by reducing it to a hxed point problem and providing con¬ 
ditions that ensure solvability. 

Let us denote by T the operator 

r = cj ■ Ve -|- £■ , (5-2) 


where C is given by 


We assume that 


Crj = —A^t] -|- h'{0, x)r] . 
h{0, x) = 0 . 


Then, fl5.ip becomes 


TUo{9, x) -h G{Uo{9, x),x) = f{9, x) 


with G given by 


G{Uo{9, x), x) = h{Uo{9,x),x) — h'{0, x)Uo{9, x) 


(5.3) 


(5.4) 

(5.5) 


We notice that T is a diagonal operator in the Fourier basis, since it is separated in the 
sum of two parts, one of which acts only on 9 and the other acting only on x. If we 
assume that C satishes HI and H2 and we denote its eigenvalues by A„, then 

Thus, we notice that T is invertible and we can reduce 05.41) to the following hxed point 
problem: 

Uo{9, x) = -T-i G{Uo{9, x),x) + T"! f{9, x) . (5.6) 

We dehne the operator T by 

nU) = -T-i G{U{9, x),x) + r-i f{9, x) . (5.7) 

Using properties of composition of functions, see Proposition [36] in Appendix O we can 
show that for Uq, Vq G then T satishes the inequality: 

||r(t/o)-r(Uo)||p,,-^ = ||r-i h{Uo,x)-T-^ h{Vo,x)-T-^h’{0,x)Uo + T-^h'{0,x)Vo\\p.^ 
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where the Lipschitz constant of the composition with G is bounded by a constant C 
times dq. Finally, if we choose ao small enough so that T is a contraction, we obtain a 
solution Uq in a ball of radius Dq. Thus, we have proven the following result. 

Proposition 16. Assume that f is in Apj^m for p > 0, j, m G 2N, j > d, m > £/2 
(m > i/2 + 1 for Neumann boundary conditions). Let h : B x V ^ C with B C C open 
set, and let V he either of the form D1 or D2. We assume that h is analytic in B and 
C^{V) n C{V) in X and let h(0, x) = 0. 

Consider model B with either D, N, P boundary conditions and assume that h satisfies, 
respectively, BCD, BCN in case of D, N boundary conditions (depending on the boundary 
condition considered for the equation). Then, the zeroth order term of model B given by 
equation (15.Ih admits a solution Uq contained in a ball around the origin in Apj^m of 
small enough radius oq ■ 


5.2. Solution of the zeroth order term for model B’. The zeroth order equation 
for model B’ is 

(a; ■ Ve)Uo{e, x) - KU,{9, x) = f{0, x) , (5.8) 

which can be solved under the assumptions In fact, dehning the operator T 

acting on U as in fl5.2p . but with A^. instead of £, we can write fl5.8p as 

TUo{9,x) = f{9,x) , 

which can be solved, because T is invertible (using Hl’-H2’). 

Indeed, the operator T is diagonal in the Fourier basis. Let us expand Uq as 


Uq{9,x) = 

k&lA n>0 


for some coefficients Uo,k,nj in a similar way, let 

keZd\{0} n>0 


Then, we obtain: 


Uo^k,n 


fk,'. 


{2mbj ■ k + A^) 


(5.9) 


for /c 7 ^ 0 (with G M denoting the eigenvalues of — A^,) and f)o,o,n = 0. From fl5.9p we 
see that the assumptions Hl’-H2’ and the regularity of f{9,x) imply the regularity of 
Uq. Because of Hl’-H2’, we have that |27ria; • fc + A^| > for some u > 0; using the 
character of the norms in Fourier coefficients, we obtain the desired result. 
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6. Proof of Theorems [9l [TO] and [TO] for model A 

In this section we present detailed arguments that complete the proof of Theorems |9l 
m and [TO] for the case of the dissipative wave equation fl2.1l) of model A; in Section [7] we 
will present the necessary modihcations in order to have the results for models A’, B, B’. 

We start by proving the existence of an approximate solution up to prescribed orders 
fSection l6.ll) : then, we bound the operator in fl3.4p providing estimates in a parabolic 
domain (Section I6.2p and we conclude by showing the existence of a solution of (13.71) 
through a hxed point argument (Section 16.31) . 

6.1. Existence of an approximate solution up to prescribed orders. In this sec¬ 
tion we describe the construction of the approximate solution up to a prescribed order 
M, as in the statement of Theorem [TO] 

For model A described by equation (12.ip . the hrst order term cq of the expansion of the 
response solution (see (13.ip ) satishes the semilinear second order elliptic equation (13.81) . 

To perform the formal manipulations that lead to the approximate solution up to order 
M, we hnd it convenient to write the equation (13.31) as 

[e{u-VeY + {uj-'V0)-eA^]Us{9,x)-eA^co{x)+eh{co{x) + Us{9,x),x) =ef{9,x) . (6.1) 

We assume that a solution cq for (13.81) can be found as described in Appendix [B] Next, 
we write formally 74 = Ue{9) in powers of e (see (13.21) ). We now show that we can dehne 
an approximate solution of (16.11) as a hnite truncation of (13.21) up to order M. Inserting 
(j3.2l) into (16. ip . we get 

OO 

'^e^[e{uj-Vef + {uj-Ve)-eA^]Uj{9,x)-eA^co{x)+eh{co{x)+Ue{9,x),x)-ef{9,x) = 0 . 

( 6 . 2 ) 

Hence the hrst order in £ in (16.2p is given by (13.91) . Since cq satishes (13.81) . then the 
equation for the hrst order in e becomes 

(w ■ Ve)TO(0, x) = f{9, x) - {f){x) (6.3) 

and it is easy to see that by the non-resonance condition (14.8p . equation (16.3p has a 
solution in Ap'j^m for some p' < p'- a proof of this fact can be found in |CCdlL13] for the 
case of the varactor equation (see also |CFG14j ) and can be straightforwardly extended 
to the present situation. In fact, let us dehne g{9,x) = f{9,x) — {f){x); let us expand 
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Ui and g as 


Ui{e,x) = 


^n{.x) Ui^k,n , g{9,x)= Y gk,n 


kfzzd n>0 


k£Zd\{0} n>0 


for suitable coefficients Ui^k,ni 


, 9k,n- From fl6.3p we obtain that 



9k,n 


which is well dehned thanks to fld.Sp . The appearance of the small divisors is the origin 
of the loss of analyticity domain. 

Note that Ui{9,x), as a solution of fl6.3p . has a free parameter, namely (17i)(a;), which 
is determined at the subsequent order. 

Recalling the dehnition fl3.5p of £, the order in fl6.2l) is 


{u ■ Ve)U2i9, x) = -[(a; • + C]Ui{9, x) , 


which admits a solution U2{9,x), provided that the average of the right hand side is 
zero. In particular, the average of Ui must satisfy the equation C{{Ui)) = 0 and since C 
satishes H1-H2, then (Ui) = 0. 

Now, by using the non-resonance condition fl4.8p up to order M, one can proceed 
recursively to compute the functions Uj in Apj^m up fo order M. We have that for any 
N < M — 1 the function Un+i satishes a recursive relation of the form (13.121) . Again, we 
must require that the average of the right hand side of fl3.12l) is zero, which provides the 
average oi Un, i.e. 


C{{Un)) = -{Sn{co,U^...,Un-i)) . 


(6.4) 


The existence of {Un) satisfying (16.4p is guaranteed by assumptions H1-H2, since the 
spectrum of C is bounded away from zero. 

In this way we obtain the approximate solution up to order M under the assumption 
(14.81) . It follows that we obtain a well dehned approximate solution to all orders under the 
condition (14.91) : for instance one can adapt the argument given in Appendix H of [CG12] . 
Finally, if we assume further that / is a trigonometric polynomial of degree J > 0, then 
LAr is a trigonometric polynomial of degree NJ and we can obtain the formal solution 
up to any order N. This concludes the proof of Theorem [101 

6.2. Bounds on the operator N^. In this section we will obtain bounds on in 
appropriate spaces, when e is contained in the domain VLb dehned as the union over a 
of the domains (14.3p . We will also remark that it is not possible to obtain the same 
bounds when e is on the imaginary axis: indeed, we present separate arguments that 
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lead us to conjecture that the bound on the spectrum of iVe cannot be obtained when e 
is imaginary. 

Note that, since C acts on the x variable only and a; ■ Ve on the 9 variable, we can 
apply separation of variables and obtain that the spectrum of is 

^n,k — (27rca * /u) “1“ (cj * /u) “1“ • 

As already pointed out, the invertibility of will follow from the fact that its spectrum 
is bounded away from zero. For a general operator, the bounds on the inverse would 
need not only to estimate the spectrum, but also the spectral projections, though this is 
trivial in this case since £ is self-adjoint and ca • Ve is anti self-adjoint, so that £ + cj-Vg 
is a normal operator. 

We study the spectrum of iV^, when e ranges in the domain 

^a,B,a = {s = ^ + iT] ^ C : ^ > B , a < \e\ < 2a} . (6.5) 

Afterwards, we will £x a in such a way that we can use the hxed point argument of 
Section 16.31 and it will turn out that the best choice is a = 2, thus leading to dehning 
the solution in the domain as dehned in 04.31) . To study the spectrum of 

Ag, we will use the maximum principle for Xn^k as a function of e] hence, we will get a 
lower bound on |A„^fc| on the boundary of the domain. 

Let us suppose that —?■ oo whenever n —)■ oo. We will show later that this as¬ 
sumption can be relaxed to encompass the case that the sequence of the A„’s has a hnite 
supremum, even if this case does not appear in the applications we have in mind. Thus, 
we can £x A G N large enough so that Xk — 1 > Ai; we will hrst provide a bound for 
1 < n < K on the whole region Xla,B,a- 

For every n G N, we want to estimate inffc^^dyio} ^ ^ ^a,B,a] therefore, we 

study the behavior of 

kA„^fc(e)P = \ — e{27r00 ■ ky + i{27100 ■ k) + eXnl"^ 

= I - (^ ir]){27roo ■ kf + i{27ioo • fc) + (^ + i'r])Xn\^ 

= [-(27ra; • k^ + Xnf + [-r]{2'Koo ■ k^ + (27ra; • k) + r]Xnf ■ 

For a given n, consider the function 

rn(r, Tj) = (-r^ + Xnf + + A„) + rf , 

where r G M, ^ -|- ir; G Clearly, 

F„(27ra; • k,^,r]) = \{^ + iv)>^n,k{^ + '^v)\^ ; 


( 6 . 6 ) 
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since infi-gR r„(r, 77 ) < r„(27ra; • k,^,ri), it suffices to bound from below 

inf^gR r„(r,^,?7). 

Let us start by considering the boundary namely we consider e as 

£ = Br]'^ + ir] (6.7) 

with 7] G M\{0}, i? > 0 large enough, say B > Bq for some Bq G M_|_. Clearly, for every 
n we have that 


inf inf \e\n,k{.£)\^ > inf |r„(r, 577 ", 77 )! . 

e=Bri°‘+iri fcGZ‘^\{0} rGlR 

We recall that HI and H2 imply that inf„>i |A„| = Ai > 0. 

Obtaining lower bounds of fl 6 . 6 l) for £ of the form fl6.7p is very simple: indeed it is 
the sum of two non-negative terms, which vanish at very different places. We analyze 
carefully each of the places where one of the terms vanishes. If none of the terms vanishes, 
the lower bound is clear. 

Let us dehne the three regions 

1+ = [v^^io ^'\/A^, + 10 

I- = [—\/^—10 ^\/^, —-\/A^ + 10 ^\/a^] 

and the complement of /+ U . When r is in such complement we have 

rn(r,e,77) > {Btf{-r^ + > Cl{B7^^f\l > Cl{Btf A? 

for a suitable constant Ci > 0. When r G /+ U /_ we have 

|An ~ T^l = ~ t| + '^I — fO (2 + 10 ^)'\/A^ , 


( 6 . 8 ) 

(6.9) 


so that we obtain 


rn(r,^,i7) > [h(An - r^) + r]' 


> - 2 |r| |? 7 | |An - r 


21 

-3^ 


> r^- 2 |r| |77|(2 + 10 -^) 10 -^A„ . 

Since we have that |r| < (H-10“^)\/A))^ and > (1 —10“^)^A„, we obtain for r G J+UJ_: 

rn(r, e, h) > (1 - 10 “^)"A„ - 2(1 + 10-3) (2 + 10-3) 10-3 a| I 77 I . 

Since we are considering Ai < A^ < Ai^, the following inequality holds for all r: 

Vn{T,i,r^)>Cl{Btf ( 6 . 10 ) 

for some constant C 2 > 0, provided I 77 I is sufficiently small to satisfy the condition 
Cl{Bri^f + 2(1 + 10-3) (2 + 10-3) 10-3 I 77 I < (1 - 10-3)^ Ai . 
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We now consider the remaining parts of the boundary of the region ^lcr,B,a-, starting 
from the circle \e\^ = and we begin from the case |A„ — r^| <6 for some positive <5. 
Then, for a sufficiently small and setting e = ^ + ir], we have 

rn(r, rj) = - \nf + (1 - 2r]T){T^ - A„) + A„ 

1 Ai ( 6 . 11 ) 

for C 3 > 0 , if cr and 6 are small enough with 0 < 5 < Ai and provided 

When |A„ — r^l > 6, we have -as before- that the minimum of r„(r, rj) is reached for 
|A„ — r^l = 5 and we obtain the following bound: 

r„(r, e, n) = kr(A. - + 2rr/(A„ - r') = + (r/(A„ -r^) + rf > ■ (6.12) 

Of course on the circle |£:|^ = 4(T^ we can reason in the same way, possibly changing the 
constants. 

Let us discuss now the case n > K] again we distinguish two cases. 

For r such that |A„ — r^| < 1 we have > Xk — 1; therefore for \r]\ sufficiently small 
and for some O 4 > 0 we obtain the bound: 

^n{r,^,r]) = \e\^{Xn-T‘^f+ t‘^+ 2Tr]{Xn-T^) >T^-2\T\\r]\>^T‘^ 

> l(Ax-l)>y >C'4e^ (6.13) 

provided |^| is sufficiently small, namely 

VS- 

Finally, for |A„ — r^| > 1, n > iF, we have 

rn{T,^,V)>e\Xn-T^\^>e ■ (6.14) 

Casting together the bounds (16.91) . (16.101) . (16.1 IF (I6.12F (I6.13p . (16.Mh . there exists a 
constant C 5 > 0 , depending on Ai, such that for every n, 

r„(r,e,r/) > C 5 = C,e , (6.15) 

since we are in the domain VLa,B,a ^ {i + irj E C \ |^|>i?|? 7 |"}. This concludes the 
bounds on the spectrum of W. 
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Remark 17. (z) The bounds providing the invertihility of N;, fail when e is on the imag¬ 
inary axis. 

iii) We do not assume that the spectrum is discrete; we could take the spectrum ranging 
over any set of the real line. 

(Hi) The discussion above does not depend explicitly on the boundary conditions as¬ 
sumed for the PDE. However the boundary conditions enter through the assumption H2. 
(iv) The case in which 

sup = A < oo 

n>l 

is even simpler than the previous discussion, since in this case we can reason exactly as 
we did just for the case n < K. 

To use the fixed point argument formulated in Section 16.31 we need the bound 

for some constant Cq. This bound will be used in fl6.18p below and, in view of fl6.15p . it 
amounts to requiring 

lel > ^6^^ 

for some constant Ce > 0. This inequality is in turn implied by 

> ^6^' , 

which is possible only if 

\Bv^\ >a(5V“ + ^') • 

Therefore, since |? 7 | < 1, we must have a < 2; in conclusion, we take a = 2, being the 
best possible exponent, thus leading to define the domain fla,B as in fl4.3p . 

6.3. Existence of the fixed point. As we have discussed in Section [3l we can rewrite 
(13.71) as a fixed point equation, namely 

U(e,x) = Np{f{0,x) - {f)(x)) - NpG{U){0,x) . 

where U denotes a function of e defined by 14 = Us{0,x). In this way, we define an 
operator T acting on functions analytic in e, taking values in given by 

T{U) = NpiJ - (/» - NpGiU) . (6,16) 

For a fixed e, we find a fixed point of T by considering a domain V C with 

T('P) C "P on which T is a contraction. Since we want to obtain analyticity in e, we 
reinterpret fl6.16p as an operator acting on a space of analytic functions in e and we 
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consider the domain V in the space Ap^j^rn,cT,B consisting of analytic functions of e taking 
values in with £ ranging on the domain VL„^b- 

We endow Apj^rn,a,B with the supremum norm 


— sup 


pj,m ? 


(6.17) 


for which Apj^rn,a,B is a Banach Space. Moreover, due to Proposition [30] of Appendix jA] 
if j > d and m > £/ 2 , then Apj^rn,a,B with the norm fl6.17p is a Banach Algebra. 

Notice that flO.lSp and the fact that is diagonal implies that we can estimate its 
norm in the domains ^la,B- We note that the inhmum is reached at the boundary of the 
domains, namely 


\\N, 


- 1 | 


\ pJ,m,(7,B — ^7 


< CjB 


(6.18) 


for some C 7 > 0, provided that B is sufficiently large. Using the Banach Algebra property 
of Apj^rn,a,B and the fact that DuG{0){9, x) = 0, we note that the operator Te is Lipschitz 
in a ball Ba{0) C Apj^m,a,B of radius a > 0. Indeed, using Proposition [36] we have that 
the Lipschitz constant of the composition with G is bounded by a constant times a. 
Thus, we have shown that 


||r(U) - nV)\\p,j^m,.,B < G,B-^a-^a\\U - • 

We continue as in [CCdlL13] by showing that T is a contraction in a ball centered 
around the approximate solution that gets mapped into itself. First, we notice that the 
approximate solution = U^{6,x) (see Dehnition [ 8 ]) satishes 

\\U^\\p,j,m,a,B < Csa 


for some Ug > 0 . 

We £x ckq > 0 which will be the radius of a ball around zero in Apj^rn,a,B, so that 
we will take the constants corresponding to this ball. We will refer to this ball as the 
ambient ball. 

Our next goal will be to identify balls around the approximate solutions such that the 
operator T maps them into themselves and is a contraction. The following discussion is 
very similar to what is done in [CCdlL13] . 

Consider a ball Bp{U^) of radius (3 around . We will impose several conditions on 
(3 that ensure that the ball is mapped into itself by T and that T is a contraction. 

The ball B^{U^) is contained in the ambient ball, Bao{0) G Apj^rn,a,B, provided that 

(6.19) 


Ggc + /3 < ocq ■ 
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Hence, we will assume fl6.19p to ensure we can use the constants of the operator in the 
ambient ball. 

The operator T is a contraction on Bp{U^) provided that 

+ /3)H-< 1 . (6.20) 


Moreover, we have that the approximate solution satishes the inequality 

for some constant Cg > 0, since is a solution at least to 0{e^) as in Theorem [9] or 
fT3l The ball BpiU^) is mapped into itself, whenever 

Cj{Csa + < (3 . ( 6 . 21 ) 


Notice that to fulhll fl6.19p . fl6.20p . fl6.2ip . we are allowed to choose (3. Namely, we want 
to show that for some B large enough and for all a sufficiently small, say a < a*{B), we 
can hud (3 > 0 such that the three conditions fl6.19p . (16.201) . fl6.2ip are satished. 

It is natural to choose 


/3 = 100 a 


and, then, (I6.19p . (I6.20p . (I6.2ip are implied by 


{Cs + 100)(j < Oq 

(6.22) 

CjiCs + 100)5“^ < 1 

(6.23) 

100C'7(C'8 + 100)5"^ + CgB-^a < 100 . 

(6.24) 


We see that we can choose B large enough so that fl6.23p is satished and, then, fl6.22p . 
fl6.24p are satished for a small enough. 

In conclusion, we obtain that T admits a hxed point in the domain V, provided a and 
B are suitably chosen. This hxed point will be a function analytic in 

As a corollary, the solution is locally unique, namely we have the following result. 


Corollary 18. For a fixed e E fla,B with a, B such that (I6.19p . (I6.20p . (I6.2ip are satisfied, 
let be an approximate solution. Then, for any 6 G we have that 

lim T” U^{e) = U{e) . (6.25) 

n—>-oo 

In particular, the convergence in (I6.25p is uniform for e G Tla,B with VLa,B as in (14.3p . 
since \\U^ — U\\ < Ca^ for a positive constant C > 0, which implies that the solution is 
analytic for e G TI^ b ■ 
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7. Modifications of the proof in Section [6] to models A’, B, B’ 

In this Section we consider models A’, B, B’, providing the necessary modihcations 
to the proof developed for model A. In particular, we concentrate on the extension of 
Theorem [10] to construct an approximate solution and on the bound of the eigenvalues 
in a suitable domain as it was done for model A in Section [621 The other parts of 
the proof can be extended to models A’, B, B’, trivially. The existence of the order zero 
solution is considered in Appendix [B] and Sections 15.1115.21 respectively. 

7.1. Model A’. For model A’ we look for a response solution of the form (13.ip and 
dehne the operators N^, C and G as 

N,U{e,x) = [{u-Vey + h^-\/e)A, + C]U{e,x) 

e 

CU{6,x) = —AxU{6,x) + h'{co{x),x)U{6,x) 

G{U){6,x) = h{co{x) + U{6, x), x) — h{co{x), x) — h'{co{x), x)U{6, x) . (7.1) 

The equation fl4.4p is equivalent to the equation 

N^Ue{0,x) + G{U^){9,x) = f{9,x) - {f){x) , 
while the function cq must satisfy (13. 8 p 

To construct an approximate solution, let us again write formally 74 as 74 = 

Then, after solving equation fl3.8p for cq (see Appendix [Bj), at the hrst order in e we need 
to solve 

(cu • Ve)A:,Ui{9,x) = f{9,x) - {f){x) , 

which yields the non-average part of 7/i using H2’ and the non-resonance condition on 
u. At the second order in e we obtain the equation 

(a; • Ve)A^U2{9, x) = -[(a; • Ve)' + C]U,{9, x) , 

from which we hrst deduce that the average (7/i) should be zero by imposing that the 
right hand side has zero average and noting that C is invertible. Then, we determine the 
non-average part of U 2 by solving the remaining equation. 

At the order iV > 3 we obtain the equation 

{uj-Ve)A^UN{9,x) = -[{uj-Vef+C]UN-i{9,x)+SN{co{x),Ui{9,x),...,UN-2{9,x)) (7.2) 

for a suitable function Sn, depending on cq and on the functions Uj with j < — 1; by 

imposing that the average of the right hand side of fl7.2p is zero we get 

/:((7/^_l)) = (4^(co,7/i...,7/Ar_2)) . 
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After fixing the average of U^-i^ we obtain the non-average part of Un by solving equation 

(El. 

To conclude the proof for model A’, we proceed to estimate the eigenvalues of the 
operator in a way similar to that of model A. 

Indeed, let us write e = ^ + it] with e belonging to the domain VL^j^b dehned in (14.3^ : 
denoting by the eigenvalues associated to —A^,, we have: 

= \ — e{2TTUj ■ kY + i{2'Ku ■ k)X^ + e\n\^ 

= ^‘^{{2'kuj ■ kY — Xnf + ■ kY + ( 2710 ; • k)X^ + V^nY ■ 

As for model A, we introduce an auxiliary function r„(r, .^, 77 ) to obtain bounds on the 
eigenvalues of N^: 

rn(r, v) = - >^nY + + An) + tX^Y ■ (7.3) 

Again we £x A' G Z such that Ax — 1 > Ai and consider hrst the case n < K. Dehne the 
regions /_ and /+ as in (I 6 . 8 lh In the region (/_ U I+Y we obtain 

rn{T,^,rj) > e{r^ - XnY > Cioexl > C.oexl 

for a suitable constant Ciq > 0 . 

Within the region /_ U J+ we have 

rn(r,^,ii) > (-T^r^ + rA^ + r^An)^ >r2(A^)2-2|r7| |r| lAn-r^lA^ 

> ^^(A^)^-2A^|r;| |r|(2 + 10-3)10-=^An 

> (1 - 10-=')2An(A^)2 - 2A^(1 + 10-=')A3/2(2 + lO’^) 10-=" \r]\ > , 

for some constant Cu > 0 , if 

CiiY^ + 2A^(1 + 10-^)(2 + 10-2) 10-2 X%^\t]\ < (1 - 10-2)2(Af)2 Ai . (7.4) 

In the case n > K, we recall the expression of r„(r, rj) in fl7.3l) . Then, we consider the 
subcase |A„ — r^| < 1, which provides > Xk — 1, so that one obtains for \ri\ sufficiently 
small: 

^n{T,^,r]) = |£|^(An-r^)^ + (A^)V^ + 2 r A^ r; (An-r^) 

> (A^)V-2 A^|r/| |r| 

> ((A^)'V>l(Af)UA«--l)>(Af)A, 

M < Ia^VAk - 1 . 


provided 
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In the case |An — t^| > 1, recalling (17.3p we obtain 

Concerning the bonndaries |e| = a and |e:| = 2(T, again we can reason as for model A. 

Note that the assumption H2’ is crucial in order to get the bound; if A^, is the standard 
Laplace-Beltrami operator, the invertibility of the operator is guaranteed only for D 
boundary conditions, because H2’ is violated for N, P boundary conditions. 

This concludes the discussion of the invertibility of the operator for model A’. The 
existence of a hxed point can be done in full analogy to model A; see Section 16.31 

7.2. Model B. To construct an approximate solution, let us write formally U^{6,x) = 
Uo{9,x) + Given the zeroth order solution (see flS.ip i as in Section l5Tl 

we proceed to determine the higher order terms, matching powers of the formal series 
expansion in the equation 

£^( 0 ; ■ Vefu,{e,x) + {uj- Ve)u,{e,x) - A,u,{e,x) + h{u,{e,x),x) = f{e,x) . ( 7 . 5 ) 

At hrst order in e, we get the equation 

{uj-\/0)Ui{e,x)-A,u,ie,x) + h'{Uo{e,x),x)u,ie,x) = o, ( 7 . 6 ) 

which can be used to determine Ui. Indeed, writing (17.6p as 

[(a; ■ V,) - A, + h'{Uo{e, x),x)] U,{e, x) = 0 , 

we may £x Ui = 0. 

Remark 19. Note that ifUo is contained in a ball around the origin in with small 

enough radius as in Proposition UR then the operator 

f = {uj-Ve)-A, + h'{Uo{e, x),x) (7.7) 

is invertible. Indeed, if we write it as the sum of the invertible operator fl5.2p . introduced 
in Section 15. il plus the multiplication operator T defined as 

T(f)= [h\UQ{6, x),x)- h'(0, x)]0 , 

which is small when Uo is in a small ball, we obtain the inverse of T = T + T by a 
Neumann series argument. 

At the generic order N > 2, we obtain the equation: 

f {UN{e,x)) = -{oo ■VefUN-2{0,x) + SN{Uo{e,x),...,UN-i{e,x)) , (7.8) 
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where Sn is a known function of the U/s with j < N. If the operator T is invertible, we 
can determine Un uniquely. 

Now we assume that it is possible to solve the zeroth order equation fIS.ip (some 
sufficient conditions have been presented is Section I5.ip as well as to solve the recursive 
equations (17.6p . (17.8p by taking Uq in a small ball around the origin in the Ap^^m norm. 
We proceed to study the conditions under which (17.5p can be solved. 

We start by introducing the operator 

Ag = e^{u: ■ Ve)^ + (a; ■ Ve) — A^. + /i^(0, x) . (7.9) 

If this operator is invertible, by the same argument as in Remark [19] the operator 

A, = ■ Vef + (u; ■ V,) - A,, + h\Uo{e, x),x) (7.10) 

is invertible whenever ||f/o||pj,m is sufficiently small. Now, if we write equation (17.5p as 

KMe,x) + H{u,){e,x) = f{e,x) ( 7 . 11 ) 

where we write = Uo + Ue with x) and 

H{Ue){6,x) = h{Ue{6,x),x) - h'{Uo{9,x),x)U^{9,x) , 

we are led to solve the equation 

U,{9,x) = -A:^[H{U,){9,x)-f{9,x)] . 

Let us dehne the operator T acting on a function U = U{9,x) by 

T[U]{9,x) ^ -A:^[H{U){9,x) - f{9,x)] . (7.12) 

Using Proposition [36] of Appendix O we can show that for U,V E Apj^m, then T satishes 
the inequality: 

\\r{U)-r{V)\Um = IIA;^ {H{U))-A:^ iH{V))\\p,,^ < Cao 1177-^11^,,,^ , 

since the Lipschitz constant of the composition with H is bounded by a constant times 
Do- 

As in the case of model A, to check that T maps a small enough ball around an 
approximate solution U^{9,x) into itself and it is a contraction, we need to investigate 
the domain on which A^ can be inverted with “good bounds”. 

The multiplier \n,k,£ associated to Ag is given by 

Xn^k,e = ■ /c)^ + 27liu) • /c + A„ , 
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where the eigenvalues An of £ = —Aa;+h'(f/o( 6 *, x),x) satisfy H1-H2 with the assumption 
that ||17o||pj,m is small. For a given n, we consider the function 

fn(r, e) = + iT + Xn 

for r G M. We are interested to evaluate the quantity inf |rn(r, £)|. This function can be 
easily analyzed geometrically, since the part corresponding to —+ zr is a parabola. 
The inhmum is generated by considering the minimum distance of the parabola from the 
quantity An (which is a real number). 

For e = 0 the parabola coincides with the vertical axis, so that if An 7 ^ 0, the distance 
is always positive. Indeed, the parabola —+ ir passes through the origin and it is 
tangent there to iM.. The axis of the parabola coincides with —e'^. 

We assume that 

Re(-£2) > 5 > 0 , (7.13) 

setting e = Brj^ + ir], then (I7.13p amounts to requiring that 77 ^ — > 0, which is 

satished for r] sufficiently small. Finally, we obtain the estimate: 


|Re(-£^r^ + ir + An)| > 5r^ + An > An , 

which ensures that the spectrum of is away from zero due to H1-H2. Therefore, 
we infer that the operator Fn(r, e) is invertible and we get uniform bounds within the 
domain 

= {e = ^ + i ?7 : Re(—e^) > 5} 

for h > 0. 

Let us conclude by considering the case of £ real which is not covered by (I7.13p . say 
e = ^ with ^ G M as in fl4.6p with 6 small enough. Then, setting 

Fn(r, 0 = I - + Anp = (An - er^ + , 


it follows that 

/Fn(r,0 = 2r(2eV2-2Ane' + l) . 

dr 

We have two cases: 

case 1. 2An^^ < 1 , so that the minimum is attained at r = 0 and one has Fn(0,^) = 

> A?; 

case 2. 2An^^ > 1 and hence the minimum is attained at t ± = (2An^^ ~ 1)/(2^^^) 

(Fn(r±,^) are equal for parity reasons) and one has 


rn(r±,0 


1 An 1 


e - 
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for ^ small enough. 

Summarizing, for e real we get 

rn(r, 0 > min{Ai, 1} . (7.14) 

7.3. Model B’. We write the solution as Ue{9,x) = Uo{9,x) + Ue{9,x), where 0^ = 
Yl'jLi^^Uj{9,x). The zeroth order solution has been already discussed in Section [521 
With respect to model B, the only modihcation is that here we do not need the assumption 
h{0) = 0. To determine the higher order terms Uj, we start by considering the equation 

■ V0fUe{9, x) + {u- Ve)U^{9, x) - A^Ue{9, x) + eh{Ue{9, x),x) = f{9, x) . (7.15) 

Inserting the series expansion for f4 info fl7.15p and matching the same powers of s, we 
get the equations for the functions Uj, j > 1. 

At the hrst order in e we obtain the equation: 

(u;-V0)[/i(9,x)-A,[/i(9,x) = -h([/o(9,x),x) . (7.16) 

Let A = (w ■ Ve) — Aa,; then fl7.16p can be rewritten as 

A[/i(9,x) = -h([/o(9,x),x) , (7.17) 

and note that the right hand side is a known function, once we solved the zeroth order 
equation. 

At the order N > 2 we get a recursive equation of the form 

AUn{9,x) = Sn{Uo{9,x),U,{9,x),...,Um-i{9,x)) , (7.18) 

for a function depending on the terms Uj, 0 < j < N, which are assumed to be 
determined at the previous steps. 

Both equations fl7.17p and (17.18!) can be solved, provided the operator A is boundedly 
invertible in the spaces Apj^m- This requirement is satished under the conditions Hl’- 
H2’ on the eigenvalues of —A^. appearing in A. 

After solving the zeroth order equation as well as fl7.18p up to a hnite order N, we 
consider the formulation of fl7.15p as a hxed point equation and establish the existence 
of solutions. 

Let us dehne the operator A^ as 

Ag = ■ Ve)^ + (cu • Ve) — A^, ; ('^•19) 

then, equation fl7.15p can be written as 

AM9,x) + H,{Ue){9,x) = f{9,x) , 
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where is defined as 

He{Ue){0,x) = eh{Ue{0,x),x) . 

We need to solve the fixed point eqnation 

U,{e,x) = -A:^[H,{U,){e,x) - f{ 9 ,x)] . (7.20) 

The invertibility of the operator in the space Apj^m has been already discnssed for 
model B and we conclude that condition fIT.lOp together with ensure that the 

spectrum of is bounded away from zero, if Uq is in a sufficiently small ball around 
the origin, as required in Proposition [T6l Comparing fl7.12p for model B and fl7.20lh we 
conclude that we can reason as for model B to apply the contraction mapping argument. 

8. Optimality of the results 

The domains of analyticity for response solutions established in Theorem [5] are not 
optimal. Clearly, many details of the argument can be optimized and it is quite possible 
that one can use better fixed point theorems or better arguments. 

Nevertheless, we want to argue in this section that the results presented cannot be 
improved very dramatically and are qualitatively optimal. 

We will present rigorous results fTheorem 1^ and heuristic arguments (Conjectures[25l 
[26l [?rp that indicate that the results obtained are qualitatively optimal and quantitatively 
almost optimal. 

In particular, we believe that the domains of analyticity of the response solutions for 
models B, B’ do not contain sectors with aperture bigger than 7r/2 for generic perturba¬ 
tions. Note that 7r/2 is precisely the critical aperture of the Phragmen-Lindelof theorem 
( [PLn8( ISZ65] ). which makes the function theoretic properties of the perturbative func¬ 
tions very tantalizing. Of course, this has deep consequences for the properties of the 
asymptotic expansions and how to recover the function from the computed asymptotic 
expansion (note that not even the uniqueness of the asymptotic expansion is clear for 
functions in these domains). 

The argument presented in this section is very general and it applies to many problems 
that can be reduced to a fixed point problem with parameters and which satisfy some mild 
conditions on analyticity and compactness. In particular, it applies to the treatment of 
the varactor problem carried out in [CCdlLl,^ , but we will not formulate here the precise 
results in this case, even if they have less technicalities than those used in this work. 
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The argument we present here goes by contradiction. We show rigorously (see Theo¬ 
rem [2T]) that, if there is a family of solutions Ue whose domain includes a resonance (see 
Dehnition l20|) . if we embed the problem in a two parameter family of problems, then for 
most families we do not have a two parameter family of solutions. The second rigorous 
result (Lemma |22]) strengthens a bit the previous one by showing that if we had solutions 
for all the one-parameter families in a neighborhood (in the space of analytic families), 
we could hnd an analytic two-parameter family. The conclusion of the two results above 
is that it will be very unlikely to hnd a family of problems, so that the domain of the 
response function includes a resonance. That is, if we hnd a family whose solutions in¬ 
clude a resonance, we can hnd arbitrarily small perturbations whose response solutions 
have a domain that does not include the resonance. 

By examining the argument carefully, and by proposing alternative points of view, 
we speculate - but we do not prove it rigorously - that this argument applies to all 
resonances simultaneously. This leads to Conjecture [251 

Of course, since our contradictions are obtained by constructing perturbations which 
cannot be continued, if we consider a restricted class of models, one has to wonder whether 
the perturbations can be constructed in this class. 

Similar lines of argumentation have appeared in the literature. Notably, we have been 
inspired by the use of uniform integrability in |Poi87a] to obtain insights on the problem 
of integrability. 

8.1. Statement of rigorous results on optimality. The key to the arguments in this 
section is the concept of resonance for a parameter family of solutions. 

Definition 20. Let he an analytic family of bounded operators from a Banach space 
to itself. We say that Eq is a resonant value for the family O^, whenever the operator 
has a zero eigenvalue. 

We say that the resonance is isolated if for all 0 < |e — £o| 1; we have that Oe is 

invertible. Note that the arguments presented here do not reguire that the resonance is 
isolated. 

In our applications to non-linear problems, say = 0, we will take as the linear 

operators O^, the derivatives at the solution = DiF^{U^). 

Of course for a general family of operators there are other alternatives between having 
an eigenvalue zero and being invertible (e.g., having continuous spectrum, having residual 
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spectrum, etc.). In our case, the operators have a spectrum which is the closure of the 
set of the eigenvalues. 

The important fact about resonances is that if Eq is a resonant value, the range of the 
operator has codimension at least 1. 

To prove our results it will be useful to introduce a two-parameter family, say 
so that it will be easier to compute obstructions generated by resonances. Precisely, our 
result is based on the following arguments: 

{i) we will show (see part {b) of Theorem [Mil that given a two-parameter family of 
problems, = 0, such that is resonant, we cannot expect to 

obtain solutions analytic in jj, near /i = 0 (we call this phenomenon “automatic 
analyticity”)] 

{ii) we will show that if there is a perturbative solution for every one-parameter family, 
there has to be a jointly analytic solution in two parameters (see Lemma [22]); 

{in) the consequence of these two results is that it is impossible that there is a solution 
that drives through the resonances for every one-parameter family (see part (a) 
of Theorem [21]) . 

We will work mainly with the equation (13.71) . but - as anticipated before - we extend 
it adding a parameter for the nonlinearity. In particular, we rewrite fl3.7p as 

M+ A^{Ue,i_i) = 0 , (8-1) 

where = G^{Ue^^) — / + (/) and is any smooth function of ja such that 

Go = G. 

We dehne the family of operators Gg as 

G, = DJ^,,o{Ue,o) = N, + Wo(G,,o) . (8.2) 

We will argue that if G^ has a resonance at £ = £o) if is very difficult to have a family 
Afj^ that allows us to have analytic in /i and which solves J^so,^{UeQ,^) = 0. 

To make all this precise, we endow the space of analytic families of linear operators 
with the topology of the supremum of the norm in a complex domain, so that it is a 
Banach space. We will always consider the domains in /i to be a ball around /i = 0. The 
domains in e could be either a ball around £ = 0 for the perturbative expansions or a ball 
around £ = Eq, where Eq is a resonance. When dealing with functions of two variables, 
we consider domains which are the product. 

The key to the argument is to show that if there are analytic solutions, the family 
has to satisfy constraints and that generic families violate them. Of course, if one 
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considers specific models in (18.Ih . it could in principle happen that the family automati¬ 
cally satishes the constraint. We will however show that this does not happen in general 
and that, even in specihc models for flS.ip . it is unlikely that one can make deformations 
satisfying the constraints imposed by the existence of analytic solutions. 

Theorem 21. Let Te be an analytic family of analytic operators from a Banach space 
to itself. Assume that there is a family Us defined in a domain of analyticity including 
Eq, such that J^eifJs) = 0 and that Eq is a resonant value for the DJ^siPs) family. 

Consider an arbitrary small ball B (Z C centered at Eq and define Ab the space of 
analytic families of operators defined in the ball endowed with the supremum topology. 

Then, we have the following results. 

(a) In any sufficiently small ball of Ab centered at iFso, can find a family of operators 
Ts such that the domain of analyticity of the solution of iFsPs) = 0 does not include Eq. 

(b) For restricted two-parameter families IFs,ii of the form (18.Ih . we have the same 
result. Namely, if there is an analytic family Usp satisfying IFsflPefl) = 0 and Eq is a 
resonant value for DIFspPefl), then for an open and dense set of families A^, we can 
find arbitrary small values /r, such that the family does not admit a solution Ue,fj, 
which is analytic near /r = 0. 

8.2. Proof of Theorem 1211 The hrst element in the proof of Theorem [2T] is the fol¬ 
lowing elementary lemma showing that if one has analytic solutions for all equations, 
then they have to be analytic in a second parameter. Afterwards, we will identify ob¬ 
structions for analyticity in two parameters near a resonance (this obstruction is very 
similar to Poincare’s obstructions to uniform integrability ( |Poi87al §81], see |dlL96j for 
a reexamination of |Poi87aj with modern techniques and for a converse of the results of 
|Poi87a] b 

Lemma 22. Consider a family of eguations PP) = 0, where is an analytic family of 
nonlinear operators. Endow the space of analytic operators with the supremum topology. 

Assume that for all in a neighborhood of IF,, in the space of analytic functions there 
is an analytic solution U^, which is locally unigue. Then, for every two-parameter family 
such that IF^ q = there exists a solution which is analytic in the two 

parameters for arbitrarily small values of fi. 

8.2.1. Proof of Lemma lEM Given a family of operators depending on two parameters 
we £x a, jd and consider the one parameter family dehned as = IFe,a£+i3- By the 
hypothesis, if a,/3 are small, we can hnd a solution 74, so that QePe) = 0. 
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Geometrically, if we let 13 vary, then the lines {e, ae + /3) form a foliation. For a different 
value of a, we obtain a transversal foliation. The solution is analytic when we restrict 
it to the leaves of each of the two transversal foliations. Note that we are using the 
hypothesis of local uniqueness to conclude that the solutions for two families are the 
same. 

Precisely, if we choose ai ^ 02 , we can consider a change of coordinates from (e,//) to 
/3i, 132 given by 


(8.3) 


CXlS — /i — /3i , (y.2^ ~ /4 — /32 , 


which gives 



—a.2j3i + q ; i /32 

^ ^- 

—CKi + 0.2 


By hypothesis the solution Up.^^ 13.2 of two-parameter family is analytic in /3i for (32 
hxed and in (32 for /3i hxed. This is the hypothesis of Hartogs theorem ( |Kranil INar71] ). 
so that we can conclude that the function is jointly analytic in /9i, (32 and, hence, 

it is jointly analytic in e, p. 

If the operators act on inhnite dimensional Banach spaces, we can reduce the proof to 
the classical result for complex valued function by observing that we can apply Hartogs 
theorem to where £ is a linear functional from the Banach space to the complex. It 

is also well known f |RS801 [RPS?] ! that functions that are analytic in this weak sense are 
strongly analytic. Alternatively, we could just note that the proof of Hartog’s theorem 
works for functions taking values in Banach spaces. □ 

Remark 23. It is amusing to note that Lemma allows one to improve the results of 
|Poi87aj . It immediately shows that if all systems in a neighborhood remained analytically 
integrable, any two parameter family would be uniformly integrable. Hence, the obstruc¬ 
tions to uniform integrability discovered by |Poi87a] show that we can get non-integrable 
systems in any neighborhood. Of course, even if Poincare was one of the creators of the 
theory of several complex variables, he did not know about Hartogs theorem. Under the 
extra assumption of uniform boundedness (which is not so unreasonable in the present 
case), the analogue of Hartogs theorem was presumably known. 

Our next result shows that there are obstructions to the existence of solutions ana¬ 
lytic in two variables in two-parameter families near resonances. This is an elementary 
application of power-series matching. Notice that the argument works in the generality 
of mappings into Banach spaces, since it is really a soft argument which applies in many 
other contexts. 
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One subtlety is that, if we consider the restricted class of families of operators as in 
it can, in principle, happen that the obstructions vanish for the restricted family. 
So, when we consider restricted families such as fIS.ip . we will need to verify that the 
family is general enough to be affected by the obstructions. 

Lemma 24. Consider the two-parameter family For some Eq, assume that the 

following equation holds: J^eo,o(Feo,o) = 0. If the range of q has codimension at least 
1 , then the space of families for which there is a solution is contained in a set of infinite 
codimension. 

Moreover, for the restricted families of the form (18.11) . if we can find f4,o solving 
IFefliUefi) = 0, then there exists an arbitrarily small pi, such that the family IFe,p does not 
have solutions close to f 4 , 0 ; which are analytic near /i = 0. 

Proof of Lemma 1241 

If there is a solution Ue,p of IFe,piUe,p) = 0 analytic in /r, we should have: 

DIF ;o,o(f4o,o) dpUsofl + IF :o,o(f4o,o) = 0 . (8-4) 

Clearly, if the perturbation is such that i9pJTo,o(f4o,o) is not in the range of DIFeoP-, 
then there is no possibility of hnding a solution of fl8.4p and, a fortiori, no possibility of 
hnding an analytic solution. 

Of course, the families IF^^p for which the hrst jet is in the range is a codimension one 
set of perturbations. Hence, the derived necessary conditions imply that the perturbation 
has to be in this set. 

Obviously, the necessary condition above is not the only one. Indeed, one can obtain 
even more obstructions for the existence of another branch by considering higher order 
terms. Matching terms up to order N, we obtain that 

+ {dpflFe„o{Ue„o) + Rn = 0, (8.5) 

where i? at is an expression involving only derivatives of order up to iV — 1. 

Clearly, the fact that Rn + (c?/x)^-^£o,o(f4o,o) range of DR^ofi gives another 

obstruction for the perturbations. 

If the range of q has codimension k, we claim that the set of families that matches 
the necessary conditions up to order iV is a submanifold of codimension Nk. In particular, 
the set of maps that satisfy all the obstructions is contained in a submanifold of inhnite 
codimension, which becomes a very meager set in the sense of Baire category theory. 
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The proof of the hrst claim of the lemma is very easy. The key observation is that 
the obstruction at order N (see (I8.5p ) involves that a given expression is in the range of 
This expression is very complicated in the coefficients of order up to — 1, but 
its dependence on the coefficient of order N is very simple. Hence, we obtain that for 
each of the functions satisfying the condition at order — 1, the obstruction at order 
N takes the form that the A^-th derivative should be an explicit expression over all the 
previous ones plus the range of Since the range of q has codimension /c, 

this increments by k the codimension of the solution of (18.hh . In the limit we obtain that 
the solution of (18.4p is contained in a set of inhnite codimension. 

The second claim of the lemma is obtained observing that for the families of operators 
as in (18.ip . then equation (18.4p gives restrictions to the derivatives Then, we 

want to show that the range of is in the complementary of the range of 

Given that in the restricted family one has and recalling that the 

set of has inhnite codimension, then we conclude that the range of is not in 

the range of and therefore the family does not admit solutions close to Uefi, 

which are analytic near /i = 0. □ 

Now we are in the position to hnish the proof of Theorem [2T1 Consider the family 
of analytic operators and let Eq be a resonant value. Let i? C C be a ball around Eq. We 
assume that all the perturbations of the family admit an analytic solution that goes across 
the ball B. If indeed there were solutions for all perturbations, then using Lemma 12^ the 
family should be analytic in two variables. However, near a resonance, which is contained 
in the ball H, by Lemma 1241 we obtain that there are many perturbations for which this 
is impossible. Hence, we conclude that the assumption that there were solutions for all 
perturbations analytic in the ball B is false. This provides part (a) of Theorem |2I1 

We conclude by mentioning that part (fe) of Theorem [21] is obtained from a straight¬ 
forward implementation of the second statement of Lemma |24l □ 

Note that the proof of Theorem |2T] goes by contradiction. We started by assuming 
that all the systems gave solutions that were analytic in a ball and we concluded that 
they were not, except in a set of inhnite codimension. This, of course, contradicts the 
hypothesis that for any perturbation, there are analytic solutions extending through a 
neighborhood of the resonant Eq and we conclude that there is one family which does not 
extend. 

Unfortunately, this does not allow us to conclude anything beyond the fact that there 
are families which do not admit solutions that extend through the resonance. Once 
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we conclude that the hypothesis fails, we cannot obtain any of the conclusions that we 
obtained from assuming its existence (and which we used to derive a contradiction). In 
particular, the argument does not allow us to conclude that the set of families which 
extend is inhnite codimension. The inhnite codimension statement was predicated on 
the fact that we had at least a solution. 

If, indeed, we could show that the set of functions for which a solution extends through 
a resonance is inhnite codimension, we could use Baire category theorem to show that 
there is a residual set of families for which the analyticity domain does not include any 
resonance. Even if the argument above does not allow us to conclude that rigorously, we 
formulate the following conjecture. 

Conjecture 25. For an open and dense set of families (in the topology indicated above), 
there is no solution defined in a neighborhood of any of the resonances. 

Notice also that for the equations considered in (I8.1jl . if we have a perturbative solu¬ 
tion of the equation, the resonances of the perturbed equation have to be close to the 
resonances of N^. Hence, we also have the following conjecture. 

Conjecture 26. Consider the problem in (18.Ih . For an open and dense set of non- 
linearities, the response solution has a singularity at a distance less than C\e\ from the 
resonances of N^. 

We hope that, perhaps, the argument used in Theorem [2T] can be strengthened to 
obtain Conjecture [26l There could also be other strategies to prove Theorem EH which 
are direct and not just by contradiction. A more constructive argument could possibly 
take the form of observing that, near the resonances, one small change in the model 
leads to a very large change in the response function. Hence, one could hope to pile up 
perturbations of the model in such a way that the model remains well dehned, but that 
the response function breaks down. 

To apply the above results to our models, one slightly delicate point is that the lin¬ 
earization depends on the solution U^. Arguing again by contradiction and in a non¬ 
rigorous way, we observe that we can compute the eigenvalues of N^Ue + fiG{Us) by using 
a perturbative expansion as in |Kat76] . Even if a full proof will be complicated, one can 
imagine that the eigenvalues can be continued analytically in p. Hence, the values of 
£ for which an eigenvalue vanishes will move continuously (of course, if there are some 
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non-degeneracy assumptions, which is reasonable to conjecture hold generically) and they 
will move differentiably. 

Hence, we are lead to the following conjecture. 

Conjecture 27. For a generic family, we can find a constant C, such that no ball of the 
form {£ : |e — £o| < with £o o, resonance for N^, is completely contained in the 

domain of analyticity of the response function. In other words, for each of the balls as 
above, we can find a point not in the domain of analyticity. 
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Appendix A. Some properties of 

A.l. Characterization of the norm in terms of the Fourier coefficients. Here 
we provide a norm equivalent to fld.lSp which can be expressed in terms of the Fourier 
coefficients. In this way, it is easy to study the boundedness of operators which are 
diagonal in the Fourier basis (products of complex exponentials in 9 and eigenfunctions 
of C in x). As before for two equivalent norms || • ||, || • ||', we will write || • || = || • ||^ 

Proposition 28. Let u G Apj^m have a Fourier expansion as in fl3.17p . We have: 
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B{k,p) 


1/2 


{{27iY \k\^ + iyx^\uk,n\^ + 


\ m 

K 


^fceZ'i\{0},nGN 

where for k we denote |/c| = |/ci| \kd\ and 


R(0,p) 


\Uo,r: 


B{k,p) = Y[a{kj,p) , 

Proof. For d = 1 we have that for k y Q-. 


a(j, p) = 


_ f 47r|j| tf j yo 


1 

47rp 


JlirikO 12 


d^e = 


—Airk Im 6 


'|Im(e)|<p 


d(Im(d)) = 


if j = 0 


347r|fc|p _ -47r|A;|p 


a{k,p) 
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Of course, when /c = 0, the integral is just Airp. For any /c 7 ^ 0 we have that the integral 
can be bounded as 


47r|fc|p 


a{k,p) 


<2nike\2^2Q < (j 


g47r|fc|p 

a{k,p) 


(A.l) 


with C- = 1 — e and C+ = 2. For d > 2 we can use Fubini’s theorem, applying (lA.ll) 
to each factor and obtaining the following inequalities: 


C^_ 


g4vr|fc|p 

B{k,p) 


< 


^27rik-9 \ 2 j2d 


< Cl 




g47r|fc|p 

B{k,p) ■ 


Finally we note that the Laplacian is diagonal in the exponentials, so that 
|(Ae + iync^ik-e^2 ^ \k\l + ly\‘^ , 


where \k \2 denotes the Euclidean norm. 

Since the exponentials are orthogonal with respect to the inner product, we obtain 

" e47r|fc|p 


Ijd 




k£l,d-\{0} 


B{k,p) 


{i 2 nr\k\i + iy\\uy 


HV. 


+ 


HV. 


5(0, p) 


This concludes the proof. □ 

Remark 29. Since the Euclidean norm \k \2 in d dimensions is equivalent to the C- 
norm \k\, we can substitute \k \2 for \k\ in the polynomial factor and obtain an equivalent 
norm in the Banach space. On the other hand, if we change the C-norm of \k\ in the 
argument of the exponential by the C-norm, we obtain a non-equivalent norm in the space 
of analytic functions. 


The space Apj^m is a closed subspace of the Sobolev space of maps from the 2d di¬ 
mensional manifold into the Banach algebra Hy. Indeed for m > i/2 (i.e., in the 
case considered in this paper), we can (i) apply Sobolev embedding theorem, (ii) obtain 
Banach algebra properties under multiplication, (iii) apply Gagliardo-Nirenberg-Moser 
inequalities for composition. 

In particular, since is a 2d-dimensional real manifold, the Sobolev embedding the¬ 
orem implies 

whenever j > d for some constant G > 0. Hence, the convergence in Hy) implies 

the uniform convergence and it is well known that uniform limits of analytic functions 
are analytic. 

Summarizing we have the following properties. 
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Proposition 30. Let p > 0, j,m E 2M. The space Apj^m of functions analytic in 9, 
endowed with the norm given in (I3.18p is a Banach algebra under multiplication, when 
j > d, m > £/2. 

Proposition 31. If we have a linear operator Ai which is diagonal in the Fourier basis, 
say 

= Xn,k , 

for suitable coefficients \n,k, then we have 

< C'sup |A„,fc| . (A.2) 

n^k 

The bound flA.2p is immediate since the operator Ai is diagonal in the Fourier basis 
and the norm is just the sum of the Fourier coefficients. 

Remark 32. We can also dispense with the assumption that C has a discrete spectrum, 
provided we assume it is self-adjoint. 

In general, if C is self-adjoint in L\(^{fD), the spectral theorem for self-adjoint operators 
in separable Hilbert spaces f |Dun58[ IHelGQl IRS801 lvN96j ) shows that there is a positive 
measure p defined on the reals and a unitary operator V from L\(j{fD) to p), where 

M is a measure space, such that 

[VCV-^f]{x)=F{x)f{x) 

for some multiplication operator F. We recall that the role of M is to account for the 
multiplicity. Very often M is just copies of R, the number of copies being the maximum 
multiplicity. The multiplication operator, then, is just a multiplication by the coordinate 
X. Of course, the unitary mapping and the measure p depend heavily on C, even if we 
do not include it in the notation. 

The above relation means that the unitary transformation sends the operator C into a 
multiplication by F{x). This allows us to define the generalized Sobolev norm as 

11 / 11 ^= [ F{xr-\Vf{x)\^dp. 

Jr 

This is guite analogous to the characterization of Sobolev spaces in terms of the Fourier 
series developed in |Taylla| . Again, by Carding’s ineguality, these spaces are eguivalent 
to the standard Sobolev spaces in the variable x and we can define Apj^m os the spaces of 
analytic functions from to H^. 
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We can also lift the spectral theorem to Fourier series. Indeed, we can consider the 
space 

xV) = L2(T'^) (g) • 

Note that the Fourier transform IF (if appropriately normalized) is a unitary operator 
from L^(T^) to = L^(M, Hence we can consider the operator 

V = H®V ■. L^{T'^)®LlciP) L2(M X M, ® dp) . 

n 

The analogue of Proposition [971 is given by the following result. 

Proposition 33. If we have a linear operator M such that 

{VMV-^f)ik,x) = X{k,x)fik,x) , 

then, we obtain 

\\M\\p,j,m<C sup |A(/c,a;)| 

A:GZ,a:Gsupport(d/i) 

for a positive constant C. 

It is important to note that the assumptions that we used in Section ldlE. giving uniform 
bounds on the operator N^, only assumed that the spectrum lied in the real line and they 
were uniform for all A real. 

The main applications of the case of non-discrete spectrum appear when T> is an un¬ 
bounded domain ofMf or another unbounded manifold. In that case, one has to take care 
of the fact that the Sobolev embeddings are different and the global regularity theory for 
elliptic eguations may be different. 


A.2. Analytic fnnctions from a Banach space to another. Now we present some 
analyticity properties of nonlinear functions in Banach spaces (the literature on this sub¬ 
ject is very wide, see for example |HP57[ Muj86| ). In analogy with the hnite dimensional 
case, there are results which show that some weak dehnitions such as differentiability 
imply stronger ones (convergence of Taylor series around every point). In the inhnite 
dimensional cases the results are more subtle since there are different notions of differen¬ 
tiability and different notions of convergence of power series, but it is true that extremely 
weak notions turn out to be equivalent to the strongest one f |HP57l Chapter III]). 

For our purposes, we only need to apply an easy implicit function theorem and to study 
the analyticity properties of the nonlinear operator G dehned in (13.6p . The Banach spaces 
in which G acts are the spaces dehned in Section 13.5.21 

The following dehnition of analytic functions will be enough for us. 
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Definition 34. Let X and Y be complex Banach spaces. We say that f ■. Vt <Z X ^ Y 
is analytic if it is uniformly differentiable at all points of namely the derivative is 
uniformly bounded and there exists a function 7 = 7 (|2:|), with 7(|-z|) —)■ 0 as \z\ —?■ 0, 
such that the following uniform bound holds: 

\\f{x + zC) - f{x) - zDf{x)C\\ < 7(1^1) (A.3) 

for all X & Q, ( & X with |^| = 1. 


Remark 35. It is clear that Definition is a rather weak notion of differentiability. 
However, it is a remarkable fact (see |HP57] . Theorem 3.17.1) that Definition fSf 
equivalent to requiring that the function f has a Taylor expansion of the form 


IS 


f(x+i) = Y,-, affine, 

^ ^ al 


that converges uniformly for ||(^|| < M for some M > 0 (indeed, < 

00 ). 

It is also true that even weaker notions of differentiability imply flA.3p . For our pur¬ 
poses, Definition will be enough, since it allows us to apply a contraction mapping 
argument. We refer to |HP57| Chapter III], |Muj86 [Her71] for other definitions of ana- 
lyticity that turn out to be equivalent. 


To make sense of the composition, we need an analytic extension of the non-linearity 
h to the complex plane with respect to its first argnment. 

Proposition 36. Let h : B x D ^ C with B an open set in C and assume that h is 
analytic in u E B and it is C^iV) fl CiV). 

Assume that dffdf^h{u,x) are bounded in B x V for ai + a 2 < m. We denote by du 
the complex derivatives and we assume that m > 1(2. 

Let Uq G such that dist(Mo(Tp, T>),C\B) > Oq with oq > 0. 

Then, for all u in a neighborhood U of Uq in we can define the operator 

Ch[u]{6,x) = h{u{e,x),x) 

from Apj^m to itself, which is an analytic operator in the sense of Definition 
Moreover, for v G Apj^m the derivative of Ch is given by 


{DCh[u]v){6,x) = h'{u{9,x),x)v{9,x) , 
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where by h' we denote the complex derivative of h with respect to its first argument, 
namely h'{u,x) = {dih){u,x) (in the proof h" will denote the second derivative of h with 
respect to its first argument). 


Proof. The proof is rather straightforward, but it requires that we interpret some ele¬ 
mentary calculations (the Taylor theorem up to order two with remainder) in different 
levels of abstraction. 

Because of Sobolev’s embedding theorem we have that the functions u, v are bounded 
and continuous, so that, for hxed 9,x and for hxed u,v, we can think of u{9,x),v{9,x) 
as numbers and assume that \v{9,x)\ is so small that u{9,x) + sv{9,x) is in the domain 
of h for s G [0,1]. 

The fundamental theorem of calculus implies that for all 9 G T^, x E V and for all 
u eU and all v G we have: 


h{u{9,x) + v{9,x),x) = h{u{9,x),x) + / h'{u{9,x) + sv{9,x),x)v{9,x) ds 


h{u{9, x),x) + h'{u{9, x),x)v{9, x) 

pS 

h"{u{9, x) + s t v{9, x),x) v{9, xY dt ds 


(A.4) 


Jo Jo 

Now we interpret the formula flA.4ll as an equality in function spaces. 

By standard Gagliardo-Nirenberg-Moser composition estimates in Sobolev spaces ( [Tayllb 
Proposition 3.9]), we have that for some Gq, > 0 depending on the norm of u, one has: 


\D°'h{u, < Ca{\\u\^icc(jd.x,)) (1 + ■ 


It is also easy to check that if u{9,x),v{9,x) are complex differentiable in 9 and h is 
also differentiable in 9, we obtain that h''{u{9,x) -f s t v{9,x),x) is a function in Apj^m 
with uniform bounds. 

Using the Banach algebra properties, we obtain indeed the desired result, since we can 


bound the integral by a constant times ||n| 




m 


the last term of flA.4j) . 


□ 


Appendix B. Solution of the zeroth order term for models A and A’ 

For models A, A’ the zeroth order term cq must satisfy fl3.8p . The literature on the 
solution of (13.81) is very wide and the results strongly depend on the form of the non¬ 
linearity. To be concrete, we quote - among the others - a result on the existence of weak 
solutions (see Proposition [37|) . which are indeed regular solutions as noticed in Remark [381 
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below, and a result on the existence of an unbounded sequence of solutions if h is odd 
(see Proposition [39]). To fix the notation, let foi^x) = {f{9,x)). 


Proposition 37. (]Pre m, Theorem 9.7) Let V C i > 3, be a bounded open set, 
f E H~^{'D); let u, u be, respectively, a lower and an wpper § solution of fld.Sp with 
D boundary conditions and with ufx) < u{x) for almost every x E V. If the function 
h : PxM —)■ M satisfies the Caratheodory conditioni\ and is increasing in its first variable. 


i.e. 


h{ui, x) < h{u 2 , x) 

with yfx) < Ui < U 2 < u{x) for almost every x E V, then fl3.8p with D boundary 
conditions has at least one weak solution u E HlifD), satisfying ufx) < u{x) < u{x) for 
almost every x eT>. 


Remark 38. For the models we are considering, it can be proved that all the weak 
solutions are smooth: under regularity conditions on the coefficients defining the elliptic 
operator A^., then a weak solution is regular. We refer the reader to Chapter f in |LU68] 
and to [AgmlOp lEvalOj . 


Assuming that h is odd, one obtains an inhnite number of solutions as shown by the 
following result. 

Proposition 39. f |Str08j . Theorem 7.2) Let T> C i > 3, be a smoothly bounded 
domain. Assume that 

i) h ; "D X M —)■ M zs continuous and odd with primitive hp{x,u) = h{x, v)dv; 

ii) there exists p <21/{I — 2) and C > 0, such that 

\h{x,u)\ < C{1 + 

almost everywhere; 

Hi) there exists q > 2, Rq > 0 such that 

0 < g hp{x,v) < h{x,v) v 

^By a lower (upper) solution of problem (13.81) . we mean a function u{u) G H^(V), such that 
h{u, ■) {h{u, •)) G {V) satisfies 

-Aa;u(x) + h{u{x),x) < (/)(x) 

{-Axu{x) + h{u{x),x) > {f){x)) 
and yfx) < 0 {u{x) > 0) when x G dV. 

function : 2? x R —>■ R satisfies the Caratheodory conditions, if h{y, •) : 2? —>■ R is measurable for 
every y G R and h{-,x) : R —>■ R is continuous for almost every x EV. 
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for almost every x, |t| > Rq; 
iv) the quantities p and q satisfy 

- 1 > ^. 

i{p — 2) Q — ^ 

Then, for any fo G the equation (13.81) with Cq{x) = 0 on the boundary ofV, has 

an unbounded sequence of solutions c^k £ k eN. 


Remark 40. i) In the above results, the Laplacian can typically be replaced by any second 
order uniformly elliptic operator with smooth coefficients (compare with |AM07j ). 

ii) The multiplicity of solutions can be studied using Lusternik-Schnirelman theory, 
which allows one to find critical points of the variational functional on a given manifold. 
The number of critical points is related to the Lusternik-Schnirelman category of the 
manifold, for which a lower bound is provided by the cup length of the manifold ( [AMOT] ). 

Hi) In our discussion we considered D boundary conditions; results are available also 
for N boundary conditions (see, e.g., |TW05[ ITan02l IWT06] ) or can be extended to P 
boundary conditions. 
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